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PREFACE 


This report presents a set of governing coupled differential equations 
for a model representing a Hybrid Heavy Lift Airship (HHLA) . These equations 
serve as the basis of a numerical study aimed at determining the aeroelastic 
stability and structural response characteristics of the HHLA. These results 
will be presented in a follow on report which will represent Part II of this 
s tudy . 

The research effort reported herein was carried out in the Mechanics 
and Structures Department at UCLA by Dr. C. Venkatesan and Professor P. Friedmann 
who served as the principal investigator. 

The authors want to take this opportunity to express their gratitude to 
the grant monitor Dr. H. Miura for his numerous constructive comments and 
suggestions . 
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SUMMARY 


This report presents a set of governing coupled differential equations 
for a model representing a Hybrid Heavy Lift Airship (HHLA) . The model consists 
of a bouyant envelope, multiple rotor systems, an underslung weight and thrusters, 
all attached to a flexible supporting structure. The dynamic equations are 
written for the individual blade with hub motions, for the rigid body motions 
of the whole model and also for the flexible modes of the supporting structure. 
The purpose of these equations is to serve as the basis of a numerical study 
aimed at determining the aeroelastic stability and structural response charact- 
eristics of the HHLA. 
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1 . INTRODUCTION 


Hybrid Heavy Lift Airship (HHLA) or Hybrid Heavy Lift Helicopter (HHLH) is 
useful for providing heavy lift capability whose potential applications are 
for logging, construction, coast guard surveilance and military heavy lift 
capability. These vehicles combine a buoyant envelope lift with lift and 
control forces generated by a multiple rotor system. A rough sketch of a 
HHLA configuration is shown in Fig. 1. Such a configuration is different from 
the conventional rotorcraft which have been considered in the past. It is 
well known that the aeroelastic and structural dynamic response problems are 
crucial for the safe design of a successful rotorcraft. Therefore it is 
essential to consider the basic aeroelastic and dynamic behavior of HHLA type 
vehicles so that the potential aeroelastic instability modes and structural 
dynamic features can be simulated and identified in the design process. 

It has been established that rotary-wing aeroelasticity is inherently 
nonlinear [Ref. 1] . Aeroelastic studies performed in both industry and re- 
search organization are indicative of this aspect. Thus the correct treatment 
of a wide class of problems in this field requires a consistent development 
of a mathematical model which includes geometrically nonlinear effects, due 
to the inclusion of finite slopes in the inertia, structural and aerodynamic 
operators. It is also well known [Ref. 2] that the unsteady aerodynamic en- 
vironment in rotorcraft is complicated, Accurate mathematical models, in- 
cluding the unsteady wake effects are rarely incorporated in aeroelastic 
analyses. In HHLA type vehicles these difficulties will be further compounded 
by interference buoyant lift. Therefore, it is clear that study of some basic 
aeroelastic effects in HHLA type vehicles is important for the effective design 
of such vehicles. 

Some of the typical problems that might be encountered by the unique con- 
figuration represented by HHLA type vehicles are described below. 

(a) Isolated Blade Instabilities : These instabilities are of the flap-lag, 

flap-pitch or coupled flap-lag-torsion type and can occur both in hover and 
forward flight. Even if some existing rotor systems which are expected to 
be free of these isolated blade type instabilities, a wake excited flap-pitch 
or coupled flap-lag-torsion flutter can occur at low thrust and low inflow 
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[Refs. 2, 3 and 4] . This situation could potentially be of interest for HHLA 
type vehicles when the rotors are lightly loaded and the buoyance ratio 3 is 
large [buoyance ratio g = buoyant lift/vehicle gross weight] . Furthermore, it 
is reasonable to concentrate primarily on the hover case for HHLA type vehicles 
because the forward speed of HHLA type vehicles will be low (i.e. y < 0.20). It 
was shown [Ref. 5] that forward flight is frequently stabilizing. 

(b) Coupled Rotor/Support System Instabilities : A rotor mounted on a moving 

or flexible support system can have additional instabilities when compared to 
an isolated blade. On the ground a mechanical instability can occur known as 
ground resonance [Refs. 2 and 6] and this instability is known to be sensitive 
to the flexibility and damping of the landing gear system. In flight, the 
coupled rotor/support system can experience an aeromechanical instability usual- 
ly denoted as air-resonance [Refs. 2 and 6], All these instabilities could be 
encountered in a HHLA type vehicle, because furthermore, the buoyancy effect 
and the flexibility of the supporting structure could modify these instabilities 
in an unexpected manner. 

(c) Vibration Problems : The vibration levels in helicopters have two peaks, when 

plotted as a function of advance ratio [Ref. 6], One peak occurs at relatively 
low advance ratios and the second at high advance ratios, Since, the advance 
ratio for the HHLA is low, this type of vehicle could experience considerable 
vibration levels. Thus, it is necessary to estimate the vibration levels and 

the resulting dynamic stresses to determine the fatigue life of the structure. 

To gain a fundamental understanding of aeroelastic effects which could be 
encountered on HHLA type vehicles due to their unique features (such as buoyancy, 
multiple rotors, flexible supporting structure and underslung load), a study of 
an idealized, simple model, representative of a typical HHLA vehicle, shown in 
Figure 2, was selected. This report presents a detailed derivation of the 
equations of equilibrium which cover the dynamics of this system. 
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2. AEROELASTIC MODEL OF AN HHLA 


2.1 Introduction and Assumption 

To study the basic aeroelastlc problems which could be encountered in an 
HHLA type configuration, a typical configuration shown schematically in Fig. 2 
will be considered. The essential features of the configuration are: 

(a) A flexible supporting structure with bending stiffness Ely(x) in the 
*s- z s plane, bending stiffness EI z (x) in the Y s -X g plane (Y s coordinate 
is normal to the figure), a torsional rigidity GJ(x) and a mass distri- 
bution m(x) . 

(b) Two rotor systems capable of providing lift, each having an arbitrary 
number of .blades N, are attached rigidly to the ends of the flexible 
structure. The distance between the center line of the structure to the 
hub center for the rotor systems is h 2 . 

(c) Two masses M F ^, M F2 having inertias I F1 and I F2 respectively are attached 
to the ends of the flexible structure. These masses and inertias represent 
the helicopters. The distance between the origin 0 g fixed in the supporting 
structure to the C.G. 's of the fuselages F^ and are & F1 and 

respectively. The C.G. of the supporting structure is at a distance h g 
from the origin 0 g . Furthermore it is assumed that the C.G.'s of the 
supporting structure and fuselages are on the X-axis. 

(d) A weight is attached to the structure. Its C.G. is at a distance h^ from 
the origin 0 g . This weight can move freely or it can be locked in a fixed 
position with respect to the flexible structure. 

(e) An envelope, providing the buoyant lift and drag D g acting at its center of 
pressure, is attached to the structure. The center of pressure is at a dis- 
tance h^ from the origin 0 g . The C.G. of the envelope is at a distance h^ 
from the origin 0 g . 

(f) Concentrated axial loads P^, p T2 s i mu l ate thrusters. 

Using this model, the dynamic equations of motion for the combined system 
consisting of two rotors, flexible structure, buoyant envelope and load are 

derived. The derivation requires four ingredients: blade equations with support 

motions, equations for the flexible structure connecting the rotors, equations 
representing the forces and moments introduced by the envelope and finally a 
representation of the dynamics of the load W^. 

Certain assumptions are introduced before writing the dynamic equations for 
this system, these are given below: 
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(1) The rotor blades are assumed to be rigid with equivalent root springs 
representing the flexibility of the blade. 

(2) The rotor blades are attached to the hub with an offset e from the axis 
of rotation (Hub center) . 

(3) The blade feathering axis is preconed by an angle $ . The blade has no 
torque offset, sweep or droop. 

(4) The feathering axis coincides with the elastic axis of the blade. 

(5) The blade cross-section is symmetric and has four distinct points: 
elastic center, mass center, aerodynamic center and tension center (Fig. 3). 

(6) The structural damping in the blade is assumed to be of the viscous type. 

(7) The rotor shaft is rigid. 

(8) The rotor speed is constant. 

(9) The rotor consists of three or more blades. 

(10) Two-dimensional quasi-steady aerodynamics is used to obtain the aerodynamic 
loads. There is no reverse flow and stall. The compressibility effect 

is neglected. 

(11) The C.G. of the fuselages are on the center line of the supporting structure 
such that the individual C.G.'s lie on a straight line. 

(12) The underslung mass is rigidly attached to the structure. 

(13) The elastic deformations of the supporting structure are at least one order 
of magnitude lower than that of the blade deformation. 

(14) Flexible structure is modeled by using free-free beam modes with arbitrary 
mass and stiffness distribution. 

(15) Aerodynamic forces and moments due to the envelope are modeled by using 
the model provided in Ref. 7. Quasisteady aerodynamic theory is used for 
blade aerodynamics, and aerodynamic interference between the rotor and 
the envelope is neglected. 

Based on these assumptions, the dynamic equations of motion for the model 
are derived, using force and moment equilibrium conditions at the connecting 
points as was done in Ref. 8. 

2.2 Ordering Scheme 

When deriving equations of motion for such a multi-rotor system, a large 
number of higher order terms has to be considered. Previous research has clearly 
indicated that many higher order terms can be neglected systematically by using 
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an ordering scheme [Refs. 1,8]. Warmbrodt and Friedmann [Ref. 8] and Levin 

have, in their derivation of coupled rotor/fuselage equations, assigned in a 

judicious manner, appropriate orders of magnitude for various terms encountered 

in the coupled rotor/fuselage equations. The ordering scheme employed in this 

study follows this approach. By assuming fuselage rotations of order e many 

additional terms will appear in the coupled rotor/fuselage equations. Such 

an ordering scheme was recently used in Ref. 9. In the earlier derivations, 

the fuselage was assumed to have only rigid body degrees of freedom and the 

3/2 

orders of magnitude of the corresponding perturbed quantities are 0(e ) . In 

the present case, the fuselage/supporting structure is being considered flexible 

and orders of magnitude are assigned also to the deformation of the supporting 

2 

structure. An order of magnitude of 0(e ) is assigned to the elastic deformations 
of the structure so that this effect appears in the hub motion while at the 
same time the number of terms in the equations remains manageable. 

The basis of the ordering scheme is a small dimensionaless parameter £ 
which represents typical blade slopes due to elastic deflections. It is known 
that for helicopter blades e is in the range 

0.1 < e < 0.2 


The ordering scheme is based on the assumption that 

1 + 0(e 2 ) I 1 
2 

i.e. terms of the order of 0(e ) are neglected in comparison with unity. The 
orders of magnitude for the various parameters governing this problem are given 
below. 

x^ p^abR ^ h 2 h* h 5 ^ 

COS \, sin \’ R ’ m ‘ R ’ R ’ R ’ R * R » R ’ R ° (1) 

5 K < > ’ i < > ‘ 0<I > 
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9x, 


( ) = 


= 0(1) 


9x, 


9 ok’ 6 Gk * ° (E > 


J. Levin, "Formulation of Helicopter Air Resonance Problem in Hover with Active 
Controls", M.S. Thesis, Mechanics and Structures Department, University of 
California, Los Angeles, September 1981. 
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3. COORDINATE SYSTEMS 


In the derivation of equations of motion of the HHLA model, various 
reference coordinate systems are used. The transformation relation between 
quantities referred in the various inertial, noninertial coordinate systems 
to be established before deriving the equations of motion. The relation be- 
tween two orthogonal coordinate systems with axes X^, Y^, and Xj , Y \ , Z_. 

with and e . , £ . , e . as unit vectors along the respective axes 

xi’ yi zi X J yJ z 3 

is 


(3.1) 


ie ) 


' e 'l 

Xi| 


1' Xj 

r yi ' 

• [ V 
3 1 


[ e • 

V zl } 

1 1 

e . 1 
( zj / 


where [T ] , the transformation matrix, can be found using the Euler angles re- 
quired to rotate the j-system so as to make it parallel to i-system. 

The S-system (Fig . 2) is an inertial system whose origin 0 g is fixed at the center line 
of the supporting structure in the unperturbed state with the Z g axis vertically upwards 
passing through the center of gravity of the envelope and X s is directed aft. 

The Sl-system is a noninertial coordinate system whose origin is also 
fixed at the same point 0 g of the supporting structure. This is a body fixed coordinate 
system which moves along with the body during perturbational motion. The S- 
system and Sl-system coincide with each other in the unperturbed state of the 
model . 

The R-system is another inertial system fixed at the center 0 of the 

n 

unperturbed hub. The directions of the axes of this system are parallel to 
that of the S-system. It should be noted that in the development of rotor blade 
equations, only one general rotor system with hub motions is considered. Con- 
sequently one set of rotor coordinate systems will be defined. These definitions 
are valid for all the rotor systems in the model. The only difference that 
will occur are the different hub motions due to the relative positions of the 
hubs with respect to the origin of the S-system. This is accounted for in the 
derivation by deriving a general expression for the motion of hub center Ojj 
due to the rigid body translation and rotation and due to the elastic defor- 
mations of the supporting structure. 
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The 1-system is a body fixed system with its origin fixed at the center 

of the hub 0^ (Fig. 4) . Prior to perturbational motion the 1-system coincides 

with the R-system. It is assumed that the 1-system and Sl-system are parallel 

systems, because as pointed out earlier that the elastic deformation slopes of the 

7/2 

supporting structure are of order 0(e ). So any small rotational motion given 

to the hub fixed 1-system, due to the elastic deformation of the supporting 
structure, is assumed negligible.' 

The perturbational translational motion at the hub center 0^ due to the 
rigid body motion and the elastic deformation of the structure is written as 




R e n + R e _ + R e D 
x xR y yR z zR 


(3.2) 


and if 0^, 0 , 0^ represent the yaw-pitch-roll rotations of the structure then 
the transformation matrix [T^ can be written as 
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i e 
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(3.3) 


Rotating 2k-system is a blade fixed coordinate system which rotates with the 
k*-* 1 blade. This 2k-system is rotated from the 1-system by the azimuth angle, 
of the k 1 ^ blade (Fig. 4) about axis. The transformation matrix is 


t T 2i] 


cost^ k 

“Sin'll 


sinijJ k 

COSljj k 


0 0 


0 

0 

1 


(3.4) 
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Rotating the 2k-system by an angle -3 (precone angle) about Y 2 ^-axis and 

translating the origin to the k blade bearing by a distance e e^^j gives 

the 3k-system (Fig. 5). The X^-axis is along the elastic axis of the 

undeformed k 1 "* 1 blade. Since 3 is of the order 0(e), sin 3 - 3 and 

P P P 

cos 3 - 1. The rotation matrix is 

P 


[ t 32 ] 


1 

0 

- 3 . 


0 

1 
0 


3 

P 

0 


(3.5) 


th 


The 4k-system (Fig. 6) is fixed in the cross-section of the k blade. Trans- 
lating 3k system an amount x^ gives the 4k system at the cross-section x^ 

of the k* 1 * 1 blade prior to elastic deformation. During elastic deformation 
of the k cn blade, i.e., flap, lag and torsion, the 4k system translates and 
rotates with the cross-section. The origin of the 4k system after the deform- 
ation is given as 

(3.6) 


( \ + V e x3k + V k S y3k + W k g Z 3k 


The rotation of the 4k system is obtained by Euler angles -3^, These 
angles represent the flap-lag-torsional rotation of the k th blade at location 
x^. The sequence of rotation is flap-lag-torsion. The transformation matrix 

IT 43 ] ^ 



1 0 0 


cos?, sin£, 0 
k k 


cos3 k 0 sinB^. 

[T 43 ] " 

0 cos^ sin<|)^ 


-sinC, cosC, 0 
k k 


0 10 


0 -sin^^ coscj)^ 


0 0 1 


-sinB^ 0 cos 3j t 


Since the angles 4^, 3^ - 0(e), the transformation matrix can be simplified 

by assuming sin 0-0 and cos 0 - 1, and be written as 


[T 43^ ~ 


A Vk 


1 


_ Wk A 


(3.7) 
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In our model, we have considered the blade as a rigid blade with root springs. 
So, the relation between the translation and rotation is 

V k = X k^k and w k = " ( " X k B k ) = x k B k (3 - 8) 


To facilitate the description of the blade element aerodynamics, the 5k 
system (Fig. 7) is defined by removing the torsional twisting of the blade 
from the 4k system which gives the rotation matrix 


l T 54 ] 


1 0 0 
0 cos(-cj)^) sin(-4>^) 

0 -sin(-(})^) cos(-^) 


when ^ is small. 


it can be written as 


[*54 1 


1 0 0 

0 1 -<f> k 

0 * k 1 


(3.9) 


(3.10) 


Summary of Coordinate Systems 

Symbol Coordinate System Unit Vectors 

S Inertial system. Origin 0 e , e , e 

J ° s xs ys zs 

fixed at the center line of 
the undistrubed supporting 
structure with Z g axis pass- 
through the C.G. of the 
envelope 


SI 


R 


1 


Noninertial body fixed. ®xsl’ ®ysl’ ®zsl 

Origin at the point 0 g on the 
supporting structure 

Inertial. Fixed at the un- e_ , e_ , e„ 

Rx Ry ’ Rz 

deformed hub center 0^. S and 
R are parallel systems 

Noninertial body fixed. Ori- a yl’ ®zl 

gin at the center of rotor 
hub Oy. SI and 1 are 
parallel systems 
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2k 


Rotates with k 1 -^ blade. Ori- 
gin at the center of the 
rotor hub Oy 



4k 


5k 


Rotates with k blade. Ori- 
gin at the k*"* 1 blade pitch 


e x3k’ e y3k’ e z3k 


bearing. Preconed, x^ axis 
coincident with blade elastic 


axis in undeformed position 


Rotates with k tl * blade. Ori- e^^, e e^^ 

gin at the elastic axis of 

the deformed blade cross- 

section at a distance x. 

k 

t'Vi 

Rotates with k blade. Ori- ® x 5k’ ^y5k’ ^z5k 
gin at the elastic axis of 
the deformed blade cross- 
section at a distance x^. 

Torsional rotation of the 
blade not included 
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4. MOTION OF THE VEHICLE 


The unconstrained vehicle has six rigid body degrees of freedom and also 
has elastic deformation of the supporting structure. Before presenting the 
equations of motion for the blade and the structure, it is necessary to es- 
tablish certain kinematical relations between the vehicle motion and the hub 
displacement, because the blade inertia and aerodynamic loads are affected by 
the hub displacement. In this section, the position vector of the origin of 
the hub centers 0^ and and the rotation vector at the rotor systems, due 
to the perturbations in rigid body translation and rotation, and elastic de- 
formation of the supporting structure, is derived. Subsequently these expres- 
sions are used in writing the blade loads. 


4.1 Kinematical Relations 

The sequence of perturbational motion of the vehicle model [Fig. 8], as- 
sumed to take place, consists first of rigid body translation of the origin 0 g of 
the supporting structure, i.e. origin of the S-system, then rigid body rotation 
in the sequence yaw-pitch-roll and finally in the perturbed position, the elastic 
deformation of the structure. Referring to Fig. 8, the perturbational trans- 
lational motion occurs along X g , Y g and Z g coordinate axes. After the rigid 
body rotation, the body fixed axes system is referred as Sl-system. The elastic 
deformations occur in the Sl-system [Fig. 9]. 

During rigid body perturbational translation, the origin 0. g is moved through 
a distance, 


R 


Os 


R e + R e +R e 
xs xs ys ys zs zs 


(4.1) 


Then the model is rotated about Z g axis through an angle 0 zg representing yaw 
motion, followed by a rotation 0 about the yawed Y g axis representing pitch. 

To represent roll, a rotation 0^ g is introduced about the yawed-pitched X g axis. 
The new position of the body axis system is SI. The transformation of unit 
vectors from Sl-system to S-system is given by 
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xs 


ys 


zs 


■ 1 

CD 

¥ 

CD 

0 +0 6 1 

, A 

/ e , \ 


zs xs ys 

ys zs xs 

1 xsl I 

e 

1 

-0+0 6 

< e , > 

zs 


xs ys zs 

j y s l / 

-0 

0 

1 

\ e . / 

L ys 

xs 

- 

y zsl ^ 


(4.2) 


Figure 9 refers to the pertubed state of the model after rigid body motion. 
It is assumed that the elastic deformations of the structure occur in this 
state. The elastic deformations are bending in X s p plane, bending in 
X s i, plane and torsion about axis. The deformations are represented 

by 


i. 

V sl 

along Y j 

2. 

W sl 

along Z , 
sl 

3. 

♦si 

along 


Position vector of the C.G. of the fuselage F 2 (is point 0^) after the de- 
formation is 


"0F2 


e , + v . | t + w 
F2 xsl si | £ f2 y sl 


sl 


&F2 Z sl 


(4.3) 


where the symbol | £ refers to the value of the appropriate displacement 
F2 

at location 

The position vector of the origin of hub in rotor system 2 (i.e. point 0^) 
after the elastic deformation is 


r 0H2 W F2 ' h 2 “sl.xln^ ®Jtsl + (v sl|(, F2 ' Vsl | < , F2 > 2 ysl 

+ <h 2 + ” s l| tF2 ' S z sl < 4 " 

Due to these deformations, there is also a rotation. The rotation along axes 
at 0„„ and 0 U _ are given by 

r Z HZ 


6 OF2 6 OH2 9 si|g, F2 e xsl W sl,x| 2 , F2 e ysl T V sl,x| ilF2 e zsl 


8 , + v 


(4.5) 
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Using equations 4.1, 4.2 and 4.4, the position vector of after the 
perturbational rigid body motions and the elastic deformation, can be 
written as (in the inertial system S) 


R 


0H2 


[R xs + (£ F2- h 2 "sl.xj^ + (v sl| -M-il- 


IaT t J V v TV V J 

ix F2 2 sl l%2 zs xs ys 


+ (ho+w.,1. )(6 +6 6 )] e 


1 o' rw 1 I > 

2 sl u F2 y 2 


s xs xs 


+ [R + (^o-h 0 w 


)6„„ + (v , | „ -h 0 4>„ , I ) 


ys F2 2 si ,x| £p 2 zs v slU F2 2 "sl| %2 


+ (h.+w .| ) (-8 +6 0 )] e 

2 si | £ f2 xs ys zs' J ys 


+ [R ZS + (£ F2 h 2 W sl,x]j, F2 )( 6 ys ) + (v sl|j, F2 ~ h 2 ^sllj^ 6 xs 


+ <V-” s iU f2 )! *„ 


(4.6) 


This relation can be more compactly written as 


R„ tio = R e + R e + R e 
0H2 x xs y ys z zs 


(4.7) 


The perturbational displacement at the hub center 0 of the rotor system 1 

Hi. 

can be obtained by replacing X. in Eq. 4.6 by -X . , 

F L F L 

The velocity at 0 _ due these perturbational motions is 
nZ 


R, 


OH 2 


fi[R - h„ w , | + h„ (0 +0 0 +0 0 )] e 

1 xs 2 sl,x|£ „ 2 ys zs xs zs xs xs 


+ fi[R + X „ 0 + v . . - h„(}> , , 

ys F2 xs sl|£ F2 2 si|s> F2 


+ h 0 (-0 +0 0 +0 0 ) ] e 

2 xs ys zs ys zs ys 


+ fi[R - X „0 + w . | ] e 

zs F2 ys sl|£ F2 zs 


(4.8) 
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where ( ) indicates derivative with respect to the nondimensional parameter 
ip, (ip = fit). The acceleration at due to these perturbational motion is 


R 


fi [R - h 0 w 


0H2 * xs 

2 " 


0 w . | + h o (0 +0 0 +20 0 ^0 8 ) e 

2 sl,x|£p2 2 ys zs xs zs xs+ zs xs xs 


+ n ®*s + l F2 0 * s + 


F2 


h 2 * sl k F2 


+ h„ (-0 +0 0 + 2§ 0 +0 0 )] e 

2 xs ys zs ys zs ys zs ys 


+ fi^[R - S,„ 0 8 + w . | ] e 

zs F2 ys sl IS, F2 zs 


(4.9) 


Equations 4.8 and 4.9 can be more compactly written as 

— • 9 • 

R 


0H2 


0H2 


x xs y ys z zs 


x xs y ys z zs 


) fi 

(4.10) 

) fi 2 

(4.11) 


The perturbational velocity and acceleration at the hub center 0^ can be 

obtained by replacing & F2 in Equations 4.8 and 4 . 9 by . 

• • 

Rigid body angular velocities are 0 zg fi about Z g axis and O^fi about yawed 
Y g -axis and 0 xg fi about yawed-pitched X g axis. The angular velocity of the model 
due to rigid body rotation alone is 

^igid 


fi(0 -6 0 ) e + fi(0 +0 0 ) e 

xs ys zs xs ys xs zs ys 


+ fi(0 - 0 0 ) e e 

zs xs ys zs 


(4.12) 


Angular velocity at 0^ and 0^ due to elastic deformation is (from equation 4.5) 


E flW.ll, e xsl - W sl,x|n e ysl + V S l,x|it „ ^sP 


(4.13) 


F2 


F2 


°F2 


Combining equations (4.12) and (4.13), using equation (4,2), the angular 
velocity at the rotor hub due to the elastic deformation and due to the rigid 
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body rotation can be obtained. The angular velocity given in the s-system 
is 


= 51(0 -0 0 ) e +51(0 +0 6 ) £ +fl(0 - 0 0 )e 

0H2 xs ys zs xs ys xs zs ys zs xs ys zs 


+ fl[ »slU p2 - “ S l,xl %2 <- e z s«„V> + 'A.xk^V^s'W 1 **. 


+ fi ^sl|£ F2 6 zs "s.l,x| + V sl ,x | IL-pr, ^ 0 xs + 0 ys 9 zs^ 


F2 


F2 


ys 


+ 51 [-$ I 0 - w . | 0 + v . ] 

sl U f 2 ys sl > x U F2 xs sl > x l%2 

which can be simplified to 


zs 


cu_ u0 = 51 [ (0 -0 0 +4 1 -. 

0H2 xs ys zs sl|£. 


) e .: 


F2 


+ (0 +09 - w , | ) e 

ys xs zs sl,x|£ F2 ys 


+ (0 -00 + v , i )e ] 

zs xs ys sl,x |£ F2 zs 


(A. 14) 


The Angular acceleration is 

"nuo = ^ 2 U0 -0 0 - 0 0 + <J> i | p ) e 

0H2 xs ys zs ys zs si 1 * _ xs 

c Z 

+ (0 + 00 +00 - w . I ) e 

ys xs zs xs zs sl,xi£ F2 ys 


+ (0 - 00 - 00 + v , ) e ] 

zs xs ys xs ys sl,x|* F2 zs 


(4.15) 


The angular velocity and acceleration at 0^ F is obtained by replacing & F2 
by - in equation (4.14) and (4.15). 

Assuming the rotations due to the elastic deformation of the supporting 
structure to be small compared to those due of the rigid body rotation, causes 
the angular velocity at the hub to be a result of rigid body rotation alone. 
Thus this expression becomes 
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5 

OH 2 


(6 -6 0)1 +(6 +00 ) e 

xs ys zs xs ys xs zs' ys 


+ (0 -0 0 ) e ] 

zs xs ys zs 


and 


(4.16) 


<M = fi 2 [ (0 - 0 0 - 0 0 ) e 

OH2 xs ys zs ys zs xs 


+ 

(0 

+06+90 ) 

e 


ys 

xs zs xs zs 

ys 

+ 

(0 

- 8 8 - 0 0 ) 

e ] 


zs 

xs ys xs ys 

zs 


(4.17) 


In the equations which follow the subscript s on the rotations will be 
deleted, since they are in essence rigid body rotations. 

Equations (4.8), (4.9), (4.16) and (4.17) will be used for deriving the inertia 
and aerodynamic loads on the rotor blades. 


F 
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. EQUATIONS OF MOTION FOR THE ROTOR 


When modelling the behavior of rotors the individual blades are con- 
sidered first, the equations of motion of the individual blade is derived by 
writing the equation of dynamic equilibrium under the action of aerodynamic, 
inertia and structural loads. These equations are derived for the general 
k tb blade. Subsequently these equations are coupled with the fuselage motion 
to provide the complete set of dynamical equations of motion for HHLA vehicle 
model [Fig. 2] . 


5.1 Blade Cross-Sectional Parameters 

In the derivation of the equations of motion of the blade, certain cross- 
sectional parameters described below, are required. The 4k-system is the 
cross-sectional coordinate system. The location of any point in the cross- 
section from the origin of the 4k-system is 

y ok ^y4k + Z ok S z4k 

The principal axes of the symmetric blade cross-section are rotated from the 
4k-system by the geometric pitch angle 8^ [Fig. 6] . 

The geometrical pitch angle is defined as 


0 Gk = V (X) + 9 ok + 0 lck COS \ + 0 isk Sin \ 


where 0g^C x ) is the pretwist 
0 Q k is collective pitch 

and Q-^sk are cyclic pitch components 

If r| ^ and £ ate the principal axes coordinates, the transformation for 


y ok and Z ok is given by 
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Differentiating equation (5.2) with respect to ip yields 


ok I 

N 

O 

7? 

> = 

6 / 

( 

Gk ) V 

ok ] 

/ y ok 


Also required are expressions which define quantities involved in performing 
integration over the blade cross sectional area. Defining: 


pdA = m ; 


■SL 


pn ok dA = mX. 


X 'JL 


^ok dA = 0 


PP ok dA = I MB3 ; ff P^ok dA = ^2 ’[[ pT1 0 k ^ok dA = ° 


from these, it follows that 


jj py ok dA = mXj cos 0 Gk ; jfjf pz^ dA = mXj sin 0 f 


fl ' Py ok dA = I MB3 C ° s20 Gk + I MB2 sin20 Gk 

fj ^ pZ ok M = I MB3 Sin20 Gk + I MB2 COs20 Gk 


fl pZ ok Z ok ^ = (I MB3 " W sin0 Gk COS0 C 


In these integrals, p is the density of the material and A T is the total 
cross-sectional area of the blade. 
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5.2 Equations of Motion for the Individual Blade 

Dynamic equations of equilibrium for a blade are obtained using the 
Newtonian approach. The equations are obtained by combining the structural 
operator with the inertial, aerodynamic and structural damping loads. Since, 
the rigid, offset hinged, spring restrained model of the blade is used in 
this study, the various distributed loads are integrated over the blade 
length and then combined together to give the equation of motion. The 
various distributed loads obtained first, are described in the following 
sections . 


5.2.1 Distributed Inertia Loads on the Blade 

The distributed inertia loads on the k 1 -^ blade are obtained by first 
determining the acceleration at a general point 'P' on the blade. The loads 
per unit volume are found from D'Alembert's principle and they are integrated 
over the cross-section to give the distributed blade loads per unit length of 
the blade. 

Acceleration at a point in the blade. 

The absolute acceleration at a point 'P ' , viewed from a translating and 
rotating coordinate system with respect to an inertial frame, is given by 


a pk “ R o + r P k + ^k x r pk + “k X r pk + “k X (a) k X r P k } 


(5.6) 


where R q is the position vector of the origin of the moving coordinate system 
with respect to the inertial system. 

r ^ is the position vector of the point 'P' in the k^ blade from the origin 
of the moving reference system and is the angular velocity of the moving 
coordinate system. 

In the present rotor blade analysis, the inertial system is the R-system 
whose origin is fixed at the undeformed hub location and the the moving reference 
frame is the 2k-system whose origin moves and rotates with the blade. 

The position vector of a point 'P' on the k*"^ 1 blade is 


p2k 


e 91 + (x. +u ) e + v 
x2k k k x3k k 


-y3k 


+ w, 


k c z3k + -^ok e y4k 


+ Z ok e z4k (5.7) 
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Transforming all the unit vectors to the 2k-system 


? P 2k = g x2k [(e + x k + V “ B p W k + y ok ( "V <}) k (e p +6 k ) + VW 


+ z ok « P -vVk )] 


+ a y2k [v k + y ok + z ok ( -*k )] 


+ a z2k tB P (x k + y + w k + y ok ( *k - ? k ( yy - yw 


+ Z ok (1 “ B p B k + Vk C k )] 


(5.8) 


Applying the ordering scheme and substituting = x^f^ and = x^C^ (based 
on rigid blade approximation), equation (5.8) becomes 

r p2k = e x2k te + X k + y ok (-? k _<f> k (B p + B k^ + Z ok (-B p _B k + yk )] 

+ a y2k [x k C k + y ok + Z ok ( ’V ] 


+a z2k [B k x k + x k B p + y ok (<t> k - ? k (B p + B k }) + z ok ] 

Taking the first and second derivative of r^^ an d applying the ordering 
scheme yields 

? P 2k = ng x2k [ ' B p B k x k + “k + y ok < -^k-\ (e p +B k ) - <t, k B k- 6 k B Gk- B p § Gk > 


(5.9) 


• • 


+ z ok< A + Vk + +A + W 1 


• • 


+ fiS y2k lx A + ?ok ' A 9 Gk =' + z ok ' A ‘ ’ckk* ! 


+ ag k2k A B k + *ok < V VV*k> - 5 A + B Gk* 


+ z ok < - S p B k + 6 p*k h + «p¥k - Vok 


+ e_. <b + e. )>] 

Gk p k 


(5.10) 
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and 


"\ 


? P 2k = g x2k [ “k - 6 p\ X k + ^k* “ ? k A (e p-^k ) ~ 2 h\ ~K\ 


- e Gk (B k + V - 20 GkV 


+ z ok < A + \ c k + 2 \ ? k + Vhc + 2e Gk c k + e Gk c k >] 

+ ^ g y2k [? k \ + y ok < " 6 Gk - 20 Gk*k ~ 9 gA > + ^ 4 Gk “ $ k >] 
+ n e z2k CPp u k + \ x k + y o k < \ ~ '^k (0 k + V ~ 2 Mk " ? k\ + ®Gk 51 


+ z ok < - 6 P \ + Wk + 2 Vk 5 k + Vk 5 k 


“ 20 Gk^k + 29 Gk i; k (0 p + 0 k ) + 26 Gk e k C k 


- 0 Gk*k + e Gk^k (B p + V - 6 Gk >] 


(5.11) 


t - Vi 

The angular velocity of the k blade is 


to = co + fie , 
k 0 R zl 

where co G is the angular velocity at the hub center 0u due to the fuselage motion 


H 


th 


and it is given by equation (4.16). The angular velocity of the k blade in 
2k system is (using equation (4.16) and noting that the R and S systems are 
parallel and inertial systems) 


co,, = fie ,, [coscj;, < 0 - 0 0 > + sinijt < 0 + 0 0 >] 

Zit xzk. xx y z k y xz 


+ fte [costy, <0 + 0 0 > + sinik < -0 + 0 0 >] 

y2k k y x z k x y z 


+ ^z2k I 1 + 0 z ] 


which can be written as 


“2k = fi[ “x *x2k + “y g y2k + (1 + “z } S z2k ] 


(5.12) 
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The angular acceleration is 

oL, = ft 2 e [costk <0 -08 -80 +0 + 0 0 > + 

2k x2k k x y z yz y xz 


sin ip, <0 +00 +00 - 0 +00 >] 

ky xz xz x yz 


+ e [cosijj. <0 +00 +00 -0 + 0 0 > + 

y2k ky xz xz x yz 


sinip, <-0 +00 +00 -0 - 00 >] 

T kx yz yz y xz 


+ 'tt e 01 [0 - 0 0 - 0 0 j 

z2k z x y x y 


which can be written as 


w 2k = Q [ "x g x2k + "y g y2k + "z S z2k ] 

In terms of dimensionless derivative in time. Op) , ( ) is replaced by 
fi 2 ( ) in equation (5.6). 

The acceleration is 


(5.13) 


i p2k = q2 &o + ? P 2k + 2 " 2 k X ? P 2k + "2k x ? P 2k + "2k x ( " 2 k x ? P 2k )] 

(5.14) 

In equation (5.14), all the quantities except the first term is relative to 
the 2k-system. The first term, i.e. the acceleration of the origin of the 
2k-system is in the inertial system, as given by in equation (4.11). This 
contribution can be transformed into components parallel to 2k-system, as 
indicated below. 


"x2k 


'y2k 


[cost!, < 

R + 

R 0 

- R 0 > 

k 

X 

y z 

z y 

sini^^< 

R - 

y 

R 0 
x z 

+ Vx > l 

[cosip^< 

R - 

y 

R 8 
x z 

+ R 0 > 

z X 

sinili. < 

: -R - 

- R 0 

+ R 0 >] 

r k 

X 

y ; 

2 z y 


+ Q e [R - R 0 + R 8 ] 

z2k z y x x y 


(5.15) 
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Equation (5.15) has been obtained after applying the ordering scheme. The 
various terms in equations (5.12) and (5.15) contain the rigid body motion 
and the elastic deformation of the supporting structure. The various other 
terms in equation (5.14) are given below. These expressions are also obtained 
after applying the ordering scheme. 


The Coriolis Term 


2 “2k X *" P 2k = 20 S x2k tu y x k B k ~ x k C k " “z X k^k 


+ y ok" “y (< »k " h % + V - Sc P k + V 


- (1 + ««X-*k 0 Gk> > 


+ "ok"" (1 + C °z ) { ~K ~ 0 Gk )>] 


+ 20 S y2k [ U k “ S p B k X k + Vk " “zVk X k - VA 

+ y ok " (-i =k " ^k A + V A^k " 9 k®Gk " P p 0 Gk " “ z’ls) 

" \ A + ®Gk )5> 


• • 


+ z ok < A + Vk + Vk + Vck - “A > 1 


+ 20 Vk ' VA - Vk + VWk + 


• • 


+ y . < -oo <f>, 8 . - oo (-c - <j>. (B + 3, ) — <t>, B, 

^ok x Y k Gk y k r k p k' T k k 


- \ 8 Gk " g p 9 Gk> > 


+ "ok " - “**k " «x 8 Gk - “y ( A + *A + *vk + AW* ] 


(5.16) 
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The Angular Acceleration Term 


2k 


X f p2k = ° \2k [ Vp X k " Wk + “ Vk 


+ y . < to d>, -60 > 

J ok y T k z 


+ Z . < 01 - 01 ) >] 

ok y z T k 


+ $1 e OI [to e + oi x, - oi 3 x, - oi x, 3, 
y2k z zk xpk xkk 


+ V". <- S k - *k «p + B k » - 


- « k - c k <6 p + B k » > 


+ z . < o) (-3 - 3, ) - to > ] 

ok z v p k x 


+ g z2k [ Vk c k - V " Vk + 


y . <oi - oi (-$, ) > 
ok x y k 


+ z ok <a) x ~ \ ( ‘ e k ~ 3 P + W > ] 


(5.17) 
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The Centripetal Acceleration Term 


5 2k X (5 2k X r p2k } = Q W " (x k + *k“z + 6 + Vz 5 


+ y 0 k < - ( - ? k -VV-y-W + Vk 




+ z ok <_( " B k -e P + - w A - “ z y + “ x 


^ z (-y e k )>] 


^ a y2k 1 “ X A - “zVk + “y B P*k + Vk B k 


' VAk 


+ y ok < - a + 2 y > 


+ Z ok < *k + V*z +fl, y + Vz >] 


+ Q g Z 2k [ “x X k + w x e + Wk + Vk C k 


+ y ok <ai y + W y“ Z ' "x^k > 


+ Z ok < “x ( A - B p + Vk - “A " “zV " “x 


- « y (a> y + \ + <f> k <u z ) >1 


(5.18) 
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Combining the various terms, i.e. equations (5.11), (5.15), (5.16), (5.17), 
(5.18) the three components of the acceleration can be written as 


i = Q, {-x - 2x, 0 ) - e - 2x^£ + costy, R + sinik R 

p k K z K. K lex lev 

x2k ^ 

+ *ok ^ + \ <B p + v + 2 “z c k + Vk + Vk - K 


+ 2 o) 4 >, + 2 w 8., + 2 d>, e„, - G. 

y k y Gk T k Gk k 


- K <W - 2 Vk - Vk - 6 Gk 6 P - 2B Gk B k - 6 Gk B k> 


+ z , [8, + 0 -<(>,?, + 20 ) g + 20 ) B + CU + <0 + d). to 

ok k p k k zk zp x y T kz 


+ 2 »k + 29 Gk + 2 “ z <*k + e Gk> - B k *Vk + 2 Vk 


+ Vk + 2B Gk^k + W’ 


or in a more compact form 


' x2k 


« 2 ra„ (c) - y ok + Sk ■»“ 


x2k 


x2k J 


(5.19) 


The three terms correspond respectively to 1) the constant part of x-component 
of the acceleration over the cross-section of the blade at a distance x k , 2) 
the part dependent on y Qk and 3) the term dependent on Z Qk respectively. This 
separation facilitates integration of those terms over the cross-section of the 
blade to yield inertia forces. It also helps in identifying the orders of 
magnitude of various terms . 

The y-component of the acceleration is 


1 p y2 k = Q { ~Vk " “z X k^k + V P X k + “y B k X k ~ “z x k^k + w z e 


+ “z X k " “x e p x k - W x & k x k + C k x k + 


cos3{3 k 


Ry + sini^k 


( -R x ) + 2u k - 2B p B k x k - 2o, x x k B k 
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+ y ok !-l -2« 2 - 25 k ) 


+ 2 ok l *k + 2 V, + "y - B p“z - V'-z - “x - 26 k + 2 V k 


+ 2 *A + 2 B gA - 2 “A - V - V } 


- a 2 U< c) + y ok 4 y> + * ok 4 Z) ) 

y2k 0lC y2k ° K P y2k 


(5.20) 


The z-component of the acceleration is 


z2k 


2 r • * • 

ft 10 ) x. + co e + & x, C + w x, L - (o e - OJ x. 

x k x y k k x k k y y k 


«• •« 


+ 2 Wk + R z + B k*k 


+ *ok 1+k - h «p +B k) - A B k - «k B k + e Gk + 2 “,^k + 


+ (!)(, +4) - DC] 

y k y x k 


+ 2 ok [ - B p B k + *pVk + 2B pVk + 6 pVk - 2 6 gA 


+ 2B gA (B p +B k> + 26 GkVk - Vk + B Gk B k < W 


4 Gk - 2 “x\ - 2 “x B Gk - 2 V- B k-*Vk-AWGk> 


-i A - i y (-6 p - 6, + Vk ) + 


+ “ (-3,, - 3 + <!>,£,) - oij - w (w + *.)]> 


k p T k^k' x y y k' 


0 2 [a„^ C ^ + y , a p ^ + z , a 
- n P z2k ok P z2k ok p z2k 


(5.21) 


Transforming these acceleration components from 2k to 3k systems yields 
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lp x3k ^ (a Px2k + a Pz2k ) 

= fi 2 [a (c) + fi a n (c) + y . (a n (y) + 3 a (y) ) 
p x2k P. p z2k ° k P x2k p p z2k 

+ z (a D (z) + $ a (z) )] 
ok P x 2k P Pz2k 


? y3k " ^ ap y2k 

= ft 2 [a_ (c) + y . a< y) + z . a (z) ) 
Py 2k y ok P y2k ok p y2k 


(5.22) 


i = Q (-& a + a_ ) 
Pz3k p Px2k Pz2k 


- o 2 r ft _ (c) , 0 (c) , t o (y) , (yh 

- U [-p a + + y i C-p a + a ) 

P Px2k p z2k ° k P Pv?v P-ot, 


P Px2k 


Pz2k 


+ z (-B a (z) + a (z) )] 
ok P p x2k P z 2k 


Before evaluating the distributed inertia forces and moments, it is worth 
noting down the relative orders of magnitude of the leading terms in various 
acceleration components . They are 

„ (c) 


p x2k 
a (C) = 

Py2k 

a (C) = 
Pz2k 


0(1) 

a (y) = 
’ Px2k 

0(e) 

, a (z) = 0(e) 
Px2k 

0(e) 

X) 

^ /-N 

K> ^ 
7? ^ 

II 

0(1) 

, a (Z) = 0(e) 

Py2k 

0(e) 

a (y) = 
’ Pz2k 

0(e) 

, a (z) = 0(e 2 ) 

Pz2k 


(5.23) 


This information is useful in neglecting higher order terms even before 
evaluating the integrals to obtain inertia forces and moments. 


Distributed Inertia Forces 

The distributed inertia force per unit length is obtained from the 
D'Alembert's principle, for the k tk blade 
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- p 5 p k dA 


(5.24) 


where p is the material density of the blade. 

Substituting the various components of the acceleration and using the cross- 
sectional parameters defined earlier, the distributed inertia forces are obtained, 
in the blade fixed 3k system. The acceleration components in the 3k system 
(from eq. (5.22)) are 

*P x5k " ° 2 1 - *k - * - 2 *A - Vz + COS \ + R y 

+ y ok l5 k - W * +k<VV - *A + 2 Vk - S z + 2 *Ak 


- 2 *A ' ♦A - 29 gA - 9 gA 


+ cosik < 4), 0 + 2<j>, e + 20_, 0 + 00 + 20 0 > 

k r k y Y k y Gk y p x p y 

+ sin^ k < -<t» k e x - 2* k e x - 20 Gk 0 x + e p e y - 2 Bp 0 x > ] 

+ Z ok l \ +e k + ~^k + 20 Z (3 k +B p } + 2 *k + 20 Gk + *k 0 z 


+ 29 x » k + 9 Gk> + +A + 2<)> k C k + ♦A + 29 Gk^k 


+ 6 Gk ? k + 0 y cos *k " S x sin,J; k ] } 


(5.25) 


Q { - x k ? k + Vk + 2u k - 2 e p\ X k - 20 z x k c k + 6 Z (x k + e) 


+ cos^ k < -0 X (e p + 0 k ) * k + R y - 2x k e k e x > 


+ sin<J/, < -0 (0 + 0, ) x, - R - 2x. 0. 0 > 

Y k y p k k x k k y 


+ y ok '- 1 - 2e z - A 1 


+ 2 ok l \ - \ + 2 Vz - V 6 k> 6 z - A + 2 »A + 2 *A 


+ A 9 Gk - 29 A - 0 Gk - 6 x "A - e y “A 1 > <5 - 26) 
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^z3k 


a (x k B k + R z + B p < x k + ' + 2 : 'k'k + 2x fcV 


+ CO S *, k < ( X k +e ) f " 6 y +20 x ) + WV 2 V + 


»* • 


+ situ)). < (x, +e) (6 +26 ) + x. C (6 -20 ) + 2x. 6 - B R > 

T k k x y' k k y x k k y p y 


+ 2«k ‘*k - Vk - - Vk + 9 Gk - B p 0 k + B p 9 z 


+ cosip. <0 +26 + r (0 - 20 ) + 2C 0 > 

k x y y x k y 


+ sin* k <e y - 26 x - h ( 6 X + 28 y ) - 2^] 


+ *ok i- B p B k - B P + »p\«k - “.y W - »,¥,-Vi 


- 2B p 9 z <+k + 6 Gk> - 29 p e Gk^k 


- 28 Gk\ + 26 gA + 28 Gk B k 0 k * 0 Gk*k + Wk 


* 2 *2 * 2 

e„. - e - e 

Gk x y 


• • I 


+ cos^ < - 26 k (* k +e ck ) - 2B y (-B k 4? k G k +> t G k +G k e Gk > 


*k <9 k + 29 y> - <“ B k + Vk> (8 y - 2 V - 2B p 9 k > 


+ 8 i"* k < - 28 » k + e Gk ) + 2e x (-e k '8t k G k -H k t k +c k e Gk ). 


- V 6 y - 2 V - <- B k + W'k + 26 y> - 2B p 9 y >)1 (5 ’ 2,) 

Substituting equations (5.25 through 5.27) in equation (5.24) the distributed 
inertia loads are evaluated. The components of these loads are 


x3k 


[x k + e + 2x k 6 z + 2x k C k - R x cos^ k - R y sini|> k J 


(5.28) 
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L y3k 


= mft t x k C k + 2x k ^ k 8 z ~ C x i, + e) 0^ - x,^,, - 2u,_ + 2x 1 ,&,.3_ + x T cos0 


k k 


k"k p I Gk 


+ cos^ k < x k e p 9 x + x k B k e x + 2x k B k e x - Ry > 


+ sin\Jj, <x.J6 + x, (3. 0 + 2x, (3. 0 + R > ] 

k k p y kky kky x 


(.5.29) 


[ z3k ' 1 '*A - E z - V *k + e + 2x k 9 z + 2x k C k > 

+ cos^< B p R x - 2z k ; k 5 x + ( V e)(5 y .28 x ) - ^(e',+28 ) > 


+ sinik <M - 2x. £, 0 - (x. +e) (0 +20 ) - x, (0 -20 ) > ] 

k p y kky k x y kky x 

(5.30) 


Distributed Inertia Moment 

The distributed moment per unit length for the k*" 11 blade is also obtained 
by D'Alembert principle by taking the integral of the vector product given 
below. 


*Ik 


= II l ~ P(y ok S y4k + z ok g z4k> x a pk ] dA 


(5.31) 


The moments are also evaluated in 3k system, and are given in component form 
below. 

q T = / / ~ P l (y i - z i ) a,, - (z , - y , 0. £. + y , (j>. ) a ] dA 

i x 3k ■'ok ok k Pz3k °k ok k k ■'ok k Py3k 


q ! y 3 k II p[(z ok y ok\^k + ^k^k^ 


Px3k 


“ ( - y ok^k " y ok*A - Z oA + z ok <f> k^k ) a Pz3k ] ^ 


q I 


z3k 


Aj 


pt(_y ok ? k “ y okW Z ok^k + Z ok <f) k ? k ) a p y3k 


‘ (y ok - z okV a p x3k ] ^ 
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Substituting for the acceleration components from equations (5.25) through 
(5.27) and making use of the integrals given in equation (5.5) the components 
of the distributed inertia moments (or torques per unit span) are evaluated. 
These expressions are given below. 


q Ix3k = ° { “i C ° S0 Gk [ - X k B k - K ~ B p (x k + 6 + 2x k ? k + 2x kV 


- VA + *k 5 A + X k^k 0 z 


+ cos« k < (x k +e) (0 y -20 x ) - x k C k C0 x +20 y ) 


“ 2 Vk 0 x + B p R x + -+k R y > 


+ sinify: - (x k +e) (e x +20 y ) - x^ (0 y -20 x ) 


- 2x. 5. 0 +PE - <j>, R > ] 
k k y p y k x 


+ mx I sin0 Gk [-x k C k + x k £ k + 2u k - 23^ -29^ 


+ 0 z (x k + e) + <J» k (3 x k + x k 3 k + R z ) 


+ cost}). < R - 0 (x. 3, + x, 3 ) - 2x. 3, 0 
T k y x k k k p k k x 


+ \ X k ( ' 9 y + 29 x ) > 


+ sin^< - R x - 0 y (x k 3 p + x k 3 k ) - 2x k 3 R 0 y 


+ \ x k (9 x + 2 V H 

+ (I MB3 COS 9 Gk + I MB2 sln °Gk ) [ ~ \ + ? A + 2C k 6 k + ? k 9 k 

; - 9 Gk + 0 P «k * V 


• • 


+ \ ( - 1 - 28 z- 2 ^k ) + 
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+ cosik < - (0 +20 + c, (0 - 20 ) + 2r 0 ) > 

k x y k y x k y 


+ < 9 y - 29 „ - 0 k «* + 20 y> - 2 W > 1 


• •• 


+ (I MB3 sI " 2 9 Gk + I MB2 ^ 0 Gk> ( +k ‘ +k + V, ' VW - 28 k 


+ 2 Vk + 2 Vk + 2 Vck 


20 z 6 k - e Gk - e x cos \ " e y sin *k ] 


+ (I MB3 - W sln0 Gk COS0 Gk t" 1 - 20 z " 2 ?k ]} 


(5.32) 


q I y 3 k = 0 {mX I COS0 Gk [ ^k (x k + e + 2x k ? k + 2x k 6 z ) - 6 k ? k X k 


- C k (x k B k + R Z + e p x k } 


+ cosi|> k < -<}) k R z - ^(x k + e) (-0 y + 20 x ) > 


+ sin^< -4» k R y - C k (x k + e)(0 x + 20 y ) > ] 


+ mxj. sin0 Gk [x k +e + 2x k ? k + 2x k 0 z - R x cos^ - R sin^ k ] 


+ (I MB3 COS 9 Gk + I MB2 Sln °Gk ) ^k^k ~ ^k “ ^k^k + ^z 


+ « p +e k > <5k-*k>**'k<+k«k +2 \ B k- 29 2 H> 


- Vk <v«k> - 29 Gk<i- 5 Gk=k 


+ 2 Wck * 2 Vk B k + Vk 


+ cosik < (0 +20 + C (0 - 20 ) + 2£, 0 ) 

k kx y k y x ky 


k 0. (0 + 20 ) 

k k x y 
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- 4 >, ( 4 >. e + 2 <j».e + 20 , 0 +00 + 20 0 ) 

k k y k y Gk y p x p y 


+ sl “*k ‘ ‘ C k »y - 26 x - 4 k «, + 26 ,> - Vx 


k B. (e - 2e ) 

k k ' y x 


»k <-*k 6 x - 2 *k 6 x - 29 Gk 9 x + B p 9 y - 2B pV " 


+ (1 MB3 Sin 0 Gk + r MB2 C ° S °Gk } ^k “ B k ~ 9 p + ^k 4 k " 29 z^ B k +9 p 0 


Vz 


• • 


" 2 \ - 26 Gk (1 + V " 29 z <*k 


^k^k " 2(f> k i; k “ ^k^k _ 0 Gk^k 


- 0 costf 1 , + 0 sinik ] 
y k x k 


(I 


HB3 - W sl “ 9 Gk COs9 Gk 1 - 4 k + 9 k - 2 \ <6 k «p> - 29 z C k + 9 z 


4 \ 6 Gk + 2 »k 6 k + Vk + 29 Gk 9 k 


+ Vk - 24 , k*k 


+ cozt k < - 2f k 6 y - 2* k e y - 26 Gk 6 y - (B p + B fc ) (6, + 26 y ) > 

+ si„* k < 2* k 5 x + 2^6, + 26 Gk 6 x - < 6p 6B k )<e y - » k )>]}(5.33) 


qi z3k " “ 1 "“I “ S0 Gk <- \ “ e ‘ 2 Vk ' 2x k 9 z + E x COSi 'k + E y si *k) 

+ mxj sln0 Gk U t U k + e + 2x k c k + 2x^3 


* B k ^V’k + V’k + Vk> 
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+ cosik < -<)>. R + B. R > 
k k x k y 


+ sin^< -^R y - B k R x >] 


+ (I MB3 C ° S 9 Gk + X MB2 Sln 9 Gk )[ ~ Sc + ^k ~ \\ ~ ®z + 2 \ 9 Gk 


- 2 *A - +A ~ 26 Gk B k “ Vk - 2 ^k 


• • 


+ cos 4\ < 41.0 + 2cj) -e + 20 o , e +06 + 20 e > 

k k y k y Gk y p x p y 


• * * » 


+ sxn V - Vx - 2 ^ 0 x - 20 Gk 0 x + 0 p 0 y - 20 p 0 x > 1 


+ (I MB3 sin20 oL- + I vn»9 cos20 cvHA < A +B v + ^~Mv +2 0 , (&i, + BJ 


Gk MB 2 


Gk^ 1 ^k p k r k v k z VH k V 


• • 


+ V, + 2 * k + 2s Gk + 26 z <» k + + Vk 


* 2 Vk + Vk + 29 Gk S k + 6 gA 2 


+ B k < *k - *k + 2 Vz * (B p + V 6 z * 2B k 


+ 2 Vk + 2 Vk + «k 9 Gk - 29 z B k - e Gk > 


- Vl'tl - *k - 2B k - 6 Gk 2 


+ co»* k < -\e - e k 8 x + \G k e x > 


+ Pin* k <* k e x - e k e y + * k c k 8 y >] 


+ (I MB3 ' W cos9 Gk sln9 Gk ( “ e k + B p + 29 z B p + 2 \ + 2< V 


+ 29 z »k + 9 Gk> + 2 Vk + 26 Gk^k 


- 2 Vk - 2i =k B k 


+ cos;}). <0 — 0 > 

T k y x k 

+ sinik < - 0 - 0 ?, >] } 

r k x y k 


(5.34) 
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5.2.2 Distributed Aerodynamic Loads 

Greenberg [Ref. 10] has derived expressions for unsteady lift and moment 
on a two dimensional airfoil executing harmonic motion in a pulsating stream 
of incompressible fluid. This derivation is an extension of Theodorsen's un- 
steady aerodynamic theory [Ref. 11]. The lift and moment expressions consist 
of two contributions. The first contribution is due to circulatory flow and 
the second one is due to noncirculatory flow. Greenberg has assumed that the 
circulatory lift and noncirculatory lift are acting in the same direction, i.e. 
normal to the resultant flow. .However other researchers using this theory have 
introduced their interpretation. For example Hodges and Ormiston [Ref. 12], 
assumed that the circulatory lift acts normal to the resultant flow and the 
noncirculatory lift acts normal to the blade chord. An examination of the 
alternative mathematical expressions for the unsteady lift indicates that 
assuming the noncirculatory part of the lift to be perpendicular to the blade 
chord is somewhat more convenient. In this study it is assumed that the 
circulatory lift acts normal to the resultant . flow and the noncirculatory 
lift acts normal to the blade chord, for mathematical convenience. 

The lift and moment expression as given by Greenberg are 

L„ = 2irp bV [V a + CAT ac(k ) e ia> v t + 

C K A o o v 

[b( j - a) S + V q B] c(kg) + h c(k h ) 

+ OV 8 c(k ) e 1 ^] (5>35) 

o v+p 

2 " * ” 

L NC = ^A 13 [ h + v $ + V(a + 8) - 3bS ] (5.36) 

where is the circulatory lift 
L Nc is noncirculatory lift 

Referring to Fig. 10 

h is the vertical displacement of the axis of rotation (positive 
downward) 

a is the constant part of the angle of attack 

8 is the time varying part of the angle of attack 

V Q is the constant part of the stream velocity 
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OV e 1 ^*" is the varying part of the stream velocity 

V = V (1 + ae 1 ^) 
o 

ba is the position of the torsion axis (axis of rotation 

of the airfoil) measured from the center of the airfoil 
section. 

The total moment due to both circulatory and noncirculatory parts is 
M = irp A b 2 [ ba ii + V ba(a + 6) - Vb ( \ - a)3 - b 2 ( + a 2 )B ] 

.+ 2irp . Vb 2 (a + i) { V a + aV ac(k ) e ia) v t : - - ■ 

A 2 o o v 

+ [ b( j - a) 8 + V q 8 ] c(kg) + h c(k h ) . v 

+ aV o 3c(k v+g ) e^v 1 } (5.37) 

where M is the pitching moment about the axis of rotation (positive nose up) . 

For low frequency oscillations of the rotor blades, the reduced frequency, 
k, is low and one can introduce the assumption that the Theodorsen's lift 
deficiency function c(k) is unity. .This is equivalent to the quasisteady as- 
sumption. Furthermore from Figs. 3 and 10 one has. 

ba = -b + x A + j = x A - (5.38) 

Substituting for 'ba' from equation (5.38) and replacing c(k) by unity, the lift 
and moment equations become 

L c = 2 ttp a bV [ h + V(a + 8) + (b - * A )B 1 (5.39) 

L NC = ttp a b 2 [ h + VB + V(a + B) - (x A - |)B ] (5.40) 

2 

M = irp A b 2 { (x A - |> [ h + V (a + B) + VB - (x A - y)B]- v \ B - | 3> 

• • 

+ 2ttp a b Vx A [ h + V(a + B) + (b - X A )3 ] (5.41) 

Equations (5.39) - (5.41) can be rewritten in a modified form by replacing the 

quantity 2 tt by the incompressible lift curve slope 'a', also replacing a + 8 

by a, which represents the total effective angle of attack, where a represents 

• • ** ** 

the constant part thus a = B ; a = B and 
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(5.42) 


• • 

= p a b V [ h + Va + (b - x^) a ] 

L NC=i P A ab2[i; + V5 + ^- < X A " T> “ ] (5.43) 

M = -j p A a b 2 { (x A - y) [ h + Va + v ct - (* A - j) 

+ a p A b V x^ [ h + Va + (b - 


- . Vb - b 

a ] 2~ a - T 


x A a] 


(5.44) 


Next the various velocity components, relative to the oscillating rotor blade 
have to be identified. Let be the free stream velocity and be the 

velocity at any point on the elastic axis of the blade due to its oscillation, 
the net air flow velocity for the blade is then 


Ak 


- V 


ECk 


(5.45) 


For a rotor blade in forward flight with constant velocity, the free stream 
velocity is 

V = V cosa e - (V sina„ + v ) e (5.46) 

Ak F R Rx F R k Rz 

where V is the forward velocity of the model 

r 

a R is the angle of forward tilt of the rotor plane 
V k is the induced velocity 

Equation (5.46) can be written in terms of nondimensional quantities y and 

A t where y is the advance ratio = Vj, cosa„/QR and A. is the inflow ratio = 
k t K k 

(V F sina R + v k )/fiR. Hence 


V Ak = ^ e Rx - X k e Rz> 


(5.47) 


This velocity can be written in terms of components along the 2k-system 


V A2k = { [cos)p k < y + A k 9 y > + sin* k < y(0 y 9 x - 0 z ) -A^ > J e x2k 


+ [cos^< y(0 y 0 x - 0 z ) -A k 0 x > + sintf> k <- y - A R 0 y > ] e y2k 
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(5.48) 


+ I“< e y + 6 z 0 x> - V S z2k> 

The velocity at any point on the elastic axis due to blade deformation is 


V„ r , = R„ + r . + <o x r 
ECk Oh pk k pk 

where Rq is the velocity of the hub center 
H 


(5.49) 


ipk is the velocity of the point 'p' on the elastic axis of the blade 
as seen in the rotating reference frame 

(0^ is the angular velocity of the rotating reference frame 


and r pk = 6 e x2k + (x k + V e x3k + w k e z3k + V k e y3k 
The various terms in equation (5.49), in 2k system, are 

V = fi{( - g p e k X k + V a x2k + x k ? k a y2k + X A a z2k } 

(Based on rigid blade assumption, i.e. = x k $ k and v^ = Xj^l and 

^ A • • • 

Ro = (He 0 , [ cost, < R + R0 - R 0 > + 

x2k k x y z z y 


(5.50) 


(5.51) 


sinik < R -R0 +R9 >] 

k y x z z y 


+ e 01 [ costk <R - R0 +R0 > + 
y2k k y x z z x 


sinijj < - R - R 0 + R 0 >] 

c x y z z y 


+ e [ R - R 0 +R0 ]} 

z2k z y x x y 


(5.52) 


and 


^ x ? p2k = fi{e x2k [ -x k?k - x k?k 0 z + cos^< 0 y (g p x k + g k x k ) > 


+ Sin h K _6 x Vk + B k X k ) " 1 


+ Vk [ X k + 6 + ^z ] 
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+ e z2k [ X k ? k ^ C0S \ 6 X + S±n \ 6 y > 

- (x k +e) < cos^ k (0 y +« x 6 z ) + sin^ k (-6^40 £)>] } 

(5.53) 


Combining equations (5.51) - (5.53) 

^EC2k = Qe x2k^ ~ *k^k ~ X k ? k 9 z + u k " S p^k X k 


+ cos\|). < R + R 6 -R0 +8 (x, B + x, B, ) > 

k x y z z y y k p k k' 


+ sinik < R - R 0 + R 6 -8 (x, B + x, B, ) > 

ky xz zx x k p k k 


• • 


+ S)e y2k { x k?k + x k + e + x k 6 z + cos<|> k < R - Vz + R^ > 


+ sinik < - R - R 0 + R 0 >} 

k x y z z y 


• • • 


+ fie { x. B, + R - R 0 + R 0 

z2k k k z y x x y 


- (x k + e) < cosij^ (0 y +0^) + sin^ (-0 x +0 y 0 z ) > 


+ x k ? k < COS \ 9 x + sin \ 9 y > * 


(5.54) 


Substituting equations (5.48) and (5.54) in equation (5.45) and applying the 
ordering scheme yields 

~^2k = n{S x2k [ ¥k (1 +V " U k +B P B k X k + cos V UR ~ V 

+ sin\< -UR9 z - ^ R0 X - R y + x k 0 R (B p + B k ) > 1 

+ %2k [ _X k “ e " X k 8 z " X k^k + COS \ < ~ VR0 z “ R y “ X k R9 x > 


+ sin^ k < -pR + R x >] 
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+ ^2k [ ~ \ g k ~ X k R “ R z + yR6 y " COS Wk 0 x “ ( V e) V 

• • 

- Sin4» k < x k c k 6 y + (x^ + e) e x > ] } (5.55) 

This velocity is again transformed into components along the 5k system where 
5k system is defined as the one whose origin is fixed in the deformed blade 
elastic axis and rotated from 4k system by removing the elastic torsional 
rotation. In the 5k system, the velocity components are 


x5k 


y5k 


z5k 


fi{ - ec k - u k ■ x k e A + (6 p +B k )(_x k R - R z + pR y 

-¥^ + cos \ <tlR -V 

+ sini|> k < - liR£ k - m z ~ R y - X k R6 x + R x C k >)(5.56) 

fi[ - *k “ e “ X k ? k " X k 6 Z + COS \ < ~ ^ R6 Z " R y - 

+ sin^ k < - yR + > ] 

n [ - x k B k - A k R +-URB y - R z - x k C k (B k + Bp) 


pR? k > 


(5.57) 


" cos \ <X k C k 9 x " (x k + e) 6 y + (3 P + 8 k )(yR ~ V > 


- sin^< x k?k 9 y + (x k + e) 0 x + (8 p + " R y ) > J 

(5.58) 

For the evaluation of the unsteady aerodynamic forces and moments, the 
various velocity terms in equations (5.42) - (5.44) have to be identified. 

Figure 11 shows that 


V = -V y5k ; h = V z5k and 5 = 0 Gk + 4> k (5.59) 

Substituting these in the lift and moment expressions, the loads per unit length 
on the k c ^ blade becomes 


L Ck " Pa abV y5k V z5k + V y5k (0 Gk + ^ (b ~ X A J (0 Gk + (5.60) 
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^NGk 2 3 ^ ^ ~ (®ri. + ^1, ^ ~ Y..tU, (0/-.1, + $1,)^ 


z5k y5k ^Gk 


y5k v Gk r k' 


- < X A - ?> < S Gk + V «' 1 


(5.61) 


\ = P A ab V y5k X A [ - V z5k + V y5k < 6 Gk + +k> " < b " V (0 Gk + 


• • 


• • 


+ 2 P A b [< X A 2 i ^ V z5k V y5k^ 0 Gk*h2 V y5k ^Gk^k^ 


- <*A - T> «Ck + *k> fi2 ) 


+ I V y5k«Gk + V “ - 5 «Gk + V “ 2 > < ! - 62 > 


The drag force is 


D k-PA ab <^> IV y5k + W 


(5.63) 


The inflow angle (f)^ is 


V 


K . -1 , z5k . 

"ik = tan { V~Z ) 

y5k 


(5.64) 


According to the assumption made previously, the circulatory lift acts normal 
to the resultant flow and the noncirculatory lift acts normal to the blade 
chord. Resolving the lift and drag forces along the 5k system (Fig. 11), 
the force components per unit length are 

P A = - p A ab V y5k [ “ V z5k + V y5k (0 Gk + V " (b “ X A> C0 Gk*V n] sin *ik 
y5k. 

- 7 P A a b 2 [ V z5k “ v y5k (0 Gk + *k )fl ' V y5k (0 Gk + V 
" (X A - Y } (0 Gk + V ^ } Sin (0 Gk + V 


" P A ab ( “T > [V y5k + V z 2 5k ] COS(), ik 


(5.65) 


P A „ = p A ab V y5k [ “ V z5k + Vk^Gk^'^^ (0 Gk^k )fl 1 COS *ik 
z5k 

+ I P A a 1,2 1 * V y5k ( ®Gk + \ )!! " Vk <<J Gk + V 


(x A 2-* ^®Gk + ^ -* COS ^Gk + ^k-* 
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P A a b ( “T > [ v y5k + V z5k ] sin<t> ik 


(5.66) 


and the torsional moment, per unit length is 


• • 


Vk ' P A * b %5k *A l -V z5k + 'y5k < 8 Ck + ‘W 51 1 

+ 2 p A ab [ ‘"'A " V [V z5k " V y5k (9 Gk +4 'k ) " V y5k Ce Gk +<1, k^ 


- <>‘A-7 ) »Gk + V “ 2l 


+ ! V y5k <®Gk + V a ‘ | 2<S Gk + V ^ 1 


(5.67) 


Assuming that angle <J>_. ^ is small, the following approximations are made. 


Sin<f) ik “ ^ik = V Z 5k /V y5k 


And also 


COS< * > ik ~ 1 


sin(0 Gk + <(> k ) = sin0 Gk + <l> k cos6 


Gk 


Cos(0 Gk + V = COS0 Gk - *k sin6 Gk 

3/2 

Note that v z5k / v y5k is 0(e), 4> ik is 0(e), 8^ is 0(e' 5 ) and ( ) is 0(e ). 

Using these approximations, the force components per unit length become 

p A y5k * »A ab V z5k l, z5k - V y5k <0 Gk + V + < b -*A> < 5 ek*V n > 

- p a b ( — ) V 2 

A a y5k 

- 1 p A a » 2 [ Ksk ~ V y5k ®Gk + V “ - Vk (0 Gk + V 


<X A ‘ 7' <e Gk + *k> ” 2lUi " e GK + ♦k COS °Gk 1 


(5.68) 


45 



P Az5k = P A 3 b V y5k 1 ' V z5k + V y5k (6 Gk + V " (b ' X A } (9 Gk + V Q 1 
+ 7 °ft a » 2 ( ',5t - V y5k < 9 Gk + V “ - %k <6 Gk + V 


‘ <*A ‘ 7> <0 Gk + V fi2] Icos6 Gk - +k “^Gk 1 


(5.69) 


Where the appropriate velocity components have to be substituted in these 
equations. These velocity components are given in equations (5 . 56) -(5 . 58) . 
The aerodynamic forces per unit length of the blade, in component form are 

c 

P A = p A a bn 2 { - [ x k 2 (1 + + 26 z ) + 2x k e ] 

y5k 

" x k ^ B k * 0 Gk + \ ~ B k + 5c ^ 9 Gk + 


- 9 k «p + V + Vck* 


+ 5c (B p + B k> < °Gk + \ ' B k >] 


- x k le k < - X k R+ « (9 Gk + *k> + “ R 6 y- R k > 


+ X k R<9 Gk + +k “ B k + 6c (S Gk + V - VW+W 


+ <-‘‘ R9 y + R z ) <8 Gk + +k- B k +t k e Gk > 


+ 5c < B p + B k )< - X k E>+B k b(6 Gk + *k> 1 


-\ R <-V + e »Gk + V + RR W 


+ (pR0 - R ) < -X, R + e 8„, + pR0 - R > 
y z k Gk y z 


- X k R » <9 Gk + V 
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cos^ k [ 2 — [Xfc ( yR?k+ yR 0 z + R y )] 

- “Sc' 1 ‘Vy - Wck + K < 0 Gk + *k - «k> + Wy 1 

- ^ ( S k < pR (B p + B k ) + 6 Gk (-pRG k - pR9 z - i ) - t k m k > 


0 . < X. R - yR 0 + R > 
y k y z 


+ W «p + e k >< ' iR<B p + B k> 


- 0 Gk < - e 0 y + (B p + e k ) (PR - R x ) 


+ <"V + V * PR(B p + B k ) > 


- \ R0 y - \ “W + V BE < B P + B k> 


+ 0 < e 6 „. + pR 0 - R > 

y Gk H y z 


+ (e z 6 Gk- ^k «p + B k» < - RE < B p + B k> > 1 


- * k R<pR (B p + B k ) + 0 Gfc (-PR? t - UR 0 z - R y ) - PK k * k > 


+ (pRB y - R z ) < PR (B p + B k > - PRC k 9 Gk > 


- X k R < pR (B + 8 k ) > 


+ [s(0 Gk + V +PR0 y - Rp)<pR(Bp + B k )> 


+ " E «p + B k> b <0 Gk + V ! 


2 c 


-simp [ 


do 


k l a [x k < yR (1 + C k + V " R X > + eyR] 

- -k 2 1 K\ - 0 GkVy - i <0 Gk + \ - B k> - 0 Gk 0 k 0 * J 
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+ COS^ k 


2 

+ COS i 


, . 2 
+ sin 


9 

-x k [e k < -p R (0 Gk + + 6 Gk R x > 

+ e x < * k R - u R e y + h > 

+ c k (Bp + B k ) < - uR (6 Gk + <t> k ) > 

• 

- ®Gk < 6 ®x + (3 P + V ( " yR9 z " V > 
+ X k R6 x ' V e8 Gk + pR8 y " V 1 

- A k R < - u R (Q Gk + 4> k > + ® Gk R x > 

+ (yR0 y - R z ) < - UR <e Gk + 4> k ) > 3 


sin^ [ - 2 -f- tu R (U R ? k + U R 6 Z + R y ) 1 

-*k 2 lz V*> 

- x k [ ~ 8 x < 2VlR (& p + V * W k 0 Gk > 

+ e y < - V R (0 Gk + V > J 
+ yR0 Gk e0 y - U R (Bp + 3 k ) < UR0 Gk + UR<)> k > 

+ 2 ^ 9 Gk \ (8 P + V 1 

|/ k U k 2 e y 2 - V l 20 y (Bp + B k > UR - uRC k 0 Gk e y 3-+u 2 R 2 (Bp + B k ) 

+ pR (B p + 8 k ) < - uRC k ( 0 Gk + \> + 0 ck ( ~ R y " yR0 2 ) > 1 


U 2 0 2 + l ^ R8 Gk (8 x + C kV * UR Vx ] 


K k 1 k x 


- U R0 Gk < e0 x + (B p + B k ) (-U R 9 Z - R y> > 


'Js. < y 2 R 2 - 2 pR R > ] } 
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1_ 

2 


2 2 

p ab B 
A 


{ - x k [ (6 k - (0 Gk + J k )> < sine Gk + -<D k cose Gk > - c k e Gk sine Gk ] 
+ sin0 Gk (yR8 y - R z ) - cos* k t - x R (0 y - QJ sinQ^ 

+ [yR8 k - yR (0 Gk + *>} < sin0 Gk + 4> k cos9 Gk > 


" ^k 0 Gk sin6 Gk^ 


- sii# k [ x k (0 x + 0 ) sin0 Gk - [yR (g p + 

+ UR (6 Gk + * k ) 1 < sln6 ck + -> k cose Gk > + UR? k e Gk sinS^J } 

(5.70) 

p A z5k * p A a bf 2(x k 2 1 ®Gk + *k ' K + S. (e Gk + V ‘ VyV + °z 6 Gk 

+ J k <6 Gk + *k-«k >+8 z e Gk ) 


+ ^ I -A k R t e(6 Gk + * k ) + UR0 - R s - * k R? k 


+ a<e Gk + < 'k - V +b < 9 Gk + V 1 ' e V 


2 * 

- cosik [ - x. 0 
k k y 


+ \ I 88 (6 p + «k> + 9 Gk ( -" Ri; k - E y - W>V - >*Vk 

+ ? k <|iR u p + V .> + e 0k <- uRC k - R y - UR0 z > 

+ f' k - S k > < - yR0 k >] + epR(Bp + B k ) + \RURC k I 

- si “«v 1 \\ 

+ x k [ - yR(0 Gk + <f> k ) + R x 0 Gk - yR (6 Gk + 4> R ) £ k 

-yR0 z e Gk -yR < 8 Gk+ <i) k -0k >+ \ 0 Gk 
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- U R<t k < e Gt + V - Ct (B p + e k ) + e 2 8 Gk >] 

+ e < - WR(6 Gk + 4> k ) > + X k RyR 

- < e (0 Gk + * k > + pR9 y - R z > 

- PR b (8 0k + i k ) ) 

. + cos^ k sinip^f x^yRO^ 

- yR<yR(6 p +B k ) + 0 Gk (-yR£ k ~R y -yR0 z )-yRi; k cJ> k > 

+ < W R < 0 Gk*V " - yR0 Gk < _yR0 z ~K >] 

+ cos 2 ^ k [ - y 2 R 2 C k (3 + 3 k )]+ sin 2 ^ k [ - x fc yR0 x 

- yR < - yR(0 Gk + d> k ) + y ft > - yR Gk l > 

+ T <V bV, -* k < g k - (0 Gk + V> 

- cos^ k < yR3 k - yR (0^ +4> k ) > 

- sinip k < - yR(B p + 3 k > - yR(0g k + * k > > > cos0 Gk 

(5.71) 

It can be seen from the expressions for the aerodynamic forces that the con- 
tribution due to (0 Gk + <J> k ) associated with equation (5.40) is absent because 
of ordering scheme. Equation (5.71) can be written as 


P A. 5k = + P f5 k C ° S0 Gk 

where represents the circulatory term and pf 2 ^ represents the non- 

A z5k A z5k 

circulatory term. 

The aerodynamic moment perunit length can be written as 


(5.72) 
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x5k 


= X A P A 


z5k 


z5k 


s 


(1) , / b, (2) b* - ,'A v , 0 2 1 ,2 

+ (X A ~ 2 } P A " « (6 r.v + * 1 ,) n ?Pi ab 


2 f A 


+ \ P A afe2 \ < e Gk + <l> k > t ~ x k " e - cos^ yR? k -sin^ k MR> 


(5.73) 


p A abn 2 x A { x k 2 (0 Gk + 4> k - e k + c k (e Gk + <j> k ) - ? k (B p + e k > 


+ 9 2 0 Gk + ^ <6 Gk +<|, k- e k >+9 z 0 Gk ] 


+ x k [-X k R + e(e Gk + 4> k ) + yR6 y - R z - A k R? k 


+ e<0 Gk + \ ' B k >+ b (6 Gk + V ] " e \ R 


- cos^. [- x, 2 0 + x- [yR(B + B,) 
k k y k p k 


+ V R y - " R6 , > “ *- Rt A 


+ «k<“ R % * V >+ 8 Ck < " UEt k - \ - “ R V 


+ » k - 6 k ) (-UK t > 1 + «|lM$ +S k > + X k RpRC k ! 


- .ln* k U k + x k i-BR(8 0k + \) 


+ R k 9 Gk - mr (e Gk + * k ) ? k - we 2 e ck 


- U R < 0 <;k + A -s k > +* k e ek 


-MS' C k l® Gk + * k > -« k «, + «*> + 6 Ak > 1 


+ e < - liR te Gk + (j) k ) > + A fc RyR 


- MR < « (0 Gk + * k ) + MR6 - R,>- u* t.» ol H k )l 
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+ costp k Sim /» k [x k tiR6 y - yR < yR(B +-3 k )+0 Gk (-pR? k -R -yR6 z )-yRC k <J> k > 

+ wuCfc < VO- <e ck + V > - wR9 Gk < - V*e 2 - i y > 1 

+ cos 2 ip k [-p 2 R 2 C k CB + B k >] + sin 2 ^ k [- x kk R9 2 

- PR ' - PR (0 Gk + * k ) + R x e Gk > - PR R 2 e Gk ] 


+ i »A ■ bV (X A - 1> { - V e k ' «Gk + V 2 


COs4) k < fiRB k - yR (0 Gk + <}> k ) > 


sin^ k < ~ M R (3 + B k ) - MR (6 Gk + 4> k )>} 


7 P A afa2 I < S Gk + V 


1 2. b 2 * # 

+ 2 P A a b 2 ^ (0 Gk + * " x k ~ e “ cos \ M R ^ k ~ sin^ k uR } 

(5.74) 


These aerodynamic force and moment expressions are transformed into components 
in 3k system because the blade dynamic equations of motion are written in 
3k system. The appropriate components are 


P A " - ‘■k P A .. - B k 2 A „ 
x3k y5k z5k 


r A 

y3k y5k 

P A = P A " P k C k P A = P A 
z3k z5k y5k z5k 


q A ' = q A 
x3k x5k 


q A = ? k q A 


y3k 


x5k 


q A & k q A 

z3k x5k 


(5.75) 

(5.76) 

(5.77) 

(5.78) 

(5.79) 

(5.80) 
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5.2.3 Distributed Structural Damping Loads 

The structural damping incorporated in this analysis is of a viscous 
equivalent type. The damping forces per unit span of the blade for flap, lag 
and torsion are respectively given as 


Flap 

»D z3k - %%% 

(5.81) 

Lead-lag 

P D yJk ‘ 

(5.82) 

Torsion 


(5.83) 


No attempt is made to eliminate these terms by considering the ordering 
schemes since these terms serve the purpose of determining the effect of 
damping on stability. 

5.2.4 Rotor Blade Equations 

In this section, the individual blade equations of motion are presented. 
For the rigid, offset hinged, spring restrained blade model used in this study, 
the distributed inertia, aerodynamic and structural damping loads are integrated 
over the length of the blade and moment equilibrium at the spring restrained 
hinge is enforced. The loads due to inertia, aerodynamic and damping, integrated 
over the blade span, are given in the following sections. Finally the blade 
equations are obtained by enforcing moment equilibrium at the root of the blade. 
Inertia Loads 

The forces and moments at the blade root due to the inertia forces are 


I3k 


J rR-e 

n 


P l3k dx k 


• R-e 


<I3k f Q ( q I3k + r p3k X P l3k ) dx k 


(5.84) 

(5.85) 


where p T ,. are the distributed inertia force on the k th blade and q__, is the 
I3k I3k 

distributed inertia moment about the elastic axis. These quantities are derived 
in previous sections of this report. 
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Recall that the position vector of any point "p" on the deformed elastic 
axis of the blade is given by 


r p3k (x k + U k ) e x3k + V’k e y3k + \\ e z 3k 


(5.86) 


and is the axial displacement primarily due to geometric shortening 




<^ 2 + \ 2 > % 


' - I \ <« k 2 + O 


(5.87) 


and 


% - - *k + W 


(5.88) 


It is assumed that the inflow is constant over the disk and the pretwist of the 

blade is zero. Hence in integrations over the blade span X and 0 . remain 

Gk 

constants. Hass per unit length of the blade is also assumed to be constant. 
The components of the inertia forces at the root of the blade in 3k system 


are 


P = mS) 2 [ + (R - e) e + — 20 + 2£ 

! x 3 k 2 2 z 2 k 


- (R-e) < R cosip, + R sinik > ] 

x r k y T k 


(5.89) 


y3k 


- ■ 4i2 1 + - 5 , - « k + 2 k B P > 


■( R -e) e0 a + 2 2 [ ( k 5 k + 8 k 8 k l + (R e)x i cos8 Gk 


+ cosil). < - ( - R ~ < 0 0 + 0. 0 + 20. 0 > - (R-e) R > 

r k 2 pxkx kx y 


+ sinik < (R ^ < 0 0 + 0, 0 + 20. 0 > + (R-e) R > ] 

k 2 p y k y ky x 


(5.90) 
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?t = mQ 2 [- g - (R-e) R - B < (1+20 +2?, ) + (R-e) e > 

i z 3j t 2 k zp 2 z k 

+ cosijt < (R-e) [g R + e(0 - 26 ) ] 

k p x y x 

- — f - 1 2 v* - % - 2 v + < 5 * + *y 1 

%* *» f 

+ sint|) < (R-e) [g R - e(9 + 20 )] 

K. p y x y 

2 

- JJk f- >Vy + «, V 2 V + 5 l®y - 2 V I>] 

(5.91) 

The components of the inertia moments, after applying the ordering scheme, are 


r + P I ~ X A P I > dx v 

x3k k k X z3k k k X y3k k 




toi-R 

2 z 


- b < to£ + e + „ + i5=sii 2 ; > 

p 3 2 3 z 3. k 

+ cos( , k < xtsii b p r x - toi! 2 ; k ^ ♦ <g y - 26 x > 
+ -toi! e <e y - 2 « x ) - to£ 5k (9 k + > 

+ .i»* k 2 ei y - to^ 2 ( k e y - to>! (e x +2e y > 


e (0 + 28 ) 

2 x y 


('r-.p's *« * 

3 ~ «k « y - 2 V > J > 


0 t , , fe)! , k + feii ^ - to)! ; Z - toi! j, 

- « k + 2 <WW + -^r- 2 ¥p 

v,* ^ »A + *=£ 2 V,- toiiR y > 
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+ sinV gp e y + ^ + «=•£ 23 k e y + Ifcsl! r x >]} 


,2, . „ r (R-e) 


+ fi { mx cos9 [ - 
J. Gic 


2 B k - «-> \ 

<«-e) 2 (R-e) 2 , (R-e) 


2 . 


‘V 2 


+ (R-e) e + 

l 2 .. 


2 




+ cosip, < ( - ~^ (0 - 20 ) + (R-e) e (0 - 20 ) 

k 2 y x y X 

- c k (5 I ♦ 2 e y) -«=§ii- 2 i k 0 x 


+ (R-e) (B p R x -+ *,R y > > 


+ simp. < - (0 + 20 ) - (R-e) e (0 + 20 ) 

k 2 x y x y 

_(R-e) 2 , /q „* \ (R-e) 2 2 * ‘ 

2 y 


2 ? k ce y - 2e x ) - 


+ (R-e) <e pRy - * k R x ) > 1 

+ . sin0 Gk ^ Cfc ♦ itSli H fe - 4fi 2 « k i k ♦ B k B k ) 
- ^2|i! 2B p 6 k - 26 A + s z 

+ (R-e) eS, + iS =f i ^ \ » p + S k ) + (R-e) ^ 

2 2 

+ cosip < (R-e) R - 0 (3 + 6, ) - (R i e) 23. 0 

k y2xpk 2 kx 

2 • 

+ 2^ \ <-% + 26 *> =• 

+ sie^ < - (R-e) R x - -2^ S y (B p + 6 t ) - 26^ 


+ -^^r- * k «* + 2 V > 1 
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• « »t 


+ (a-.) [ (I„ B3 COS 2 e Gk + I MB2 slr, 2 6 Gt , ! -} k + c,e k + 2 Sfc 6 k 


+ Vk * 9 Gk + V'k - V 


+ V 1 - 2« z - 2s k ) + e k c k 


+ cosip. < - (0 + 28 +?. '( 6 - 20 ) + 2C0 ) > 

k x y k y x k y 


M • 


+ sin(f> k < - O y - 28 x - C k (8 X + 28 y ) - 25^) > ] 
+ (i MB 3 Sln29 Ck + : ME2 “'^Gk 5 [, k * \ + 2f k\ 


• • 


6 c «p + S k> - 2f! k + 2 Vk 


+ 2 Vk + 2 ^k 6 Gk - 28 A - 0 Gk 


-0 cosip - e sinip , ] 

x k y k 


+ (I MB3 " I MB2^ sin6 Gk COS0 Gk 1 “ 20 z “ (5.92) 


-y3k 


'£ % t + -A -i x3t - <\ + V ”i z3k > 8x k 


. Pt , ^2 , ia^i! + «z|i! e + toil 2 e z + 

- “L cos» k - — ‘i- R y sin» k ) } 

-- 2 ,-^5 k -^a z 

- S < fR ~^- + fR ~ y — e + — ~ 3 ? 26 + C ~^ — 25 > 

p 3 2 3 z 3 k 

+ cos4 , k < «z|)i 6p ; x _ «=§>! 25 k e x + <S^- (5 y - ^ 
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+ e < 5 y ' 2 V - 5. <§ * + 2 V > 

+ ^ t v i y . 2 - k J y _ toll ( e x + M y) 

2 • 3 

- (R ~2 } ,e(6 x + 20 y ) - ( '- ~3 } ? k (0 y - 20 x ) > ] (5.93) 


t z3k ' /„ ^zJk + <Xfc + V Pl y3k ' ^ V’ dXt 


(R-e ) 3 V | (R~e ) 3 ! . P n \ | (R~e ) 3 0 „ p 

3 C k + 3 2(C k ? k + W + 3 2f Vp 


+ costp < (R ~f - ^ 0 0+ (R ~^ ) 0 0 + (R ~^ - 20- 0 - (R ~^ ~ R > 

k 3 p x 3 k x 3 kx 2y 

+ simp, < - R ~| - - 0 0+ (R ~y~ - 0 0 + - R ~y ) ~ 20 0 + - (R ~|^— R > 
k 3 p y 3ky 3 ky 2x 

- C k ( «S) 2 [ + «z|>! . + «=12! 2 ; z + is=|ii 2 J k 

(R ~ e) R cosi|3. - (R ~ - f- - R. sirup,.] ) 


2 x ^^k 2 y ^“ T k J 

(5.94) 

The order of magnitude of leading terms in these expressions for the loads are 
listed below for convenience 


and 


x3k 


x3k 


= 0(1) 

P T = 0(e) 

y3k 

P = 0(e) 

z3k 

= 0(e 2 ) ; 

Q t = 0(e) 

y3k 

; Q = 0(e) 

z3k 


The loads at the root due to distributed aerodynamic loads are 
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and 


A3k 


X 

X 


p A3k dx k 


<A3k - / (q A3k + r p3k X P A3k } dx k 

o r 


P A B k ' [9 Gk + *k - B k + S k (9 Gk + V -«k <B p + V 

x3k 

+ B z B Gk + VVVV + B Z V 
+ I-A k R + a(0 Gk + * k ) + UR6 y - K z - \K k 


+ '<% + *k - B k> + b(6 Gk + V 1 - «- e > e \ R 


cos*.,. t- 8 y + [ GR(6 p +B k ) + 9 ( , t (-kR Ck -R y -pRe z ) 


- >"Vk + V R <B p * B k> 


+ 6 Gk ( -" K k - R y - bR V - ^'k <+k - B k» 


+ (R-e)[eyr(0 p + 3 k > + A k RyR<; k ] ] 


3 2 

sin, k 1 K * 1-«<W + R x 6 Gk-“ R ( e Gk*»k> 0 k 


UR9 z 0 Gk " yR (0 m, + K ~ K > + K Q 


Gk T k K k 


x Gk 


“ R < 'k <9 0k + V ‘ «k <B p +B k> + 6 z 6 Gk > I 


+ (R-e) [ -epR C0 Gk + *. R ) + X ^RpiR 


- yR < e (0 Gk H k ) + yR6 y - R z > 


- yR b (0 Gk + * k ) ] ] 


rR— ^ 

+ cosik sinik [ - — yR0 

k k 2 y 
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+ (R-e) L-PR < PR (B + B k ) + 20 Gk (-uRC k -R y -pR6 z ) > 

- 2pR? k 4> k ] 1 

+ cos 2 ^ [ - (R-e) y 2 R 2 £ k (B + B fc ) ] 

2 • 

+ sin 2 ^ [ - yR6 x - (R-e)yR<-yR(0 Gk +(J) k ) + 2R x 0 Gk >]} 


• • 


1 ‘V“’ 2 « 2 «k [ - - (9 Gk + 


- cos^k [ (R-e) < yR^ - yR ( 0 Qk + (J^) > ] 

- siriij^ [ (R-e) < - pR(3 + S^) - y R (0 Qk + (j^) > ] ] cos0 ( 


Gk 


, ^2 r C do r (R-e) 3 , 
(>A abn 'k [ 'T [ 3 1 


1 \ ' 9 Gk + \ - \ ) 


(R-e) 


T~ [ \ (-yR) + yt (0 Gk + * k - y 


+ ^ R6 y + V e Gk ] 


(R-e) ‘ 


+ (R-e) [ -A R R < - A r R > ] 

(R-e) 3 * 

- cos^ [ - e„o r 


3 y Gk 

T" [ ^ PR (B p + y + C-4> k + 3 k ) pR (6 p + Bk)] 
+ (R-e> [" A k RyR (B + B^ - A k RyR (B p + (3^) ] ] 


sirup 


2c do (R-e) 2 


k 1 a 
2 


pR 


I' W' 1% < (I WV > > 
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+ (R-e) [ - X k R <-yR (® Gk +\))-UR0 Gk (PRB y - R z )]] 

2 

+ cos4) k sityj^ [ - — e y ( - y R 6 Gk ) 

+ (R-e) [ - yR (g p + < y R 6 Gk + y R< f' k > ] ] 

+ cos 2 ^ k [ (R-e) y 2 R 2 ( 3 p + 3 k ) 2 ] 

2 

+ sin\ [ [ - yR 0 Gk 0 x ] 

- -f 2 ( R ~e) [ yV ] ] } 

+ | p A abV- ?k { (R-e) [ - yR0 Gk sine^ ] (- cos^)} 

(5.95) 

P A = P a abfl 2 { - -|2. [ (R ~^ (1 + 2^ + 20 ) + 2 e ] 

y3 k “ “ " A z 

- «k % + V * ®z 6 Gk > 

+ <6 p + V <e Gk + K. - 6k") 

2 # 
-■ ffi f-I% < -\ I( + e(e Gk + V + >*VV 

+ V s ' 9 Gk + \ i + ^ <0 Gk + V " \ (6 p +B k ) 

+ Vck > 

+ (- m y + V < e Gk + (v - \ * ^0 Gk > 

+ «p- V K“ «GK + VI 
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+ (R-e) [ - A R<- A R + e(6 +4>, ) + PR0 - R > 

k k Gk k y z 


+ (liRB - R ) < - A R + e 6^, + yR0 - R > 
y z k Gk y z 


- \ “> < 6 Gk + V 1 

2c , /. v2 • 

- cost}) [ — [ A - S/ ( yR£ + yR0 + R ) ] 

k a 2 k z y 

- ^ [ -Vy - C kV G k + K ‘V +V B k> + Wy> 

(R-e)" 


— [B k < PE(B p + V + 


e < A, R - yR0 + R > 
y k y z 


+ V e P + <V 4 “ E < e P + \> =• 


- V ' - »% + ,3 P + s k )( ™ - v > 


+ <- *k + V ‘ m <6 p + \ ) 3 

- \ EB y - k < e Gk + *k> «< B p +B k> 


+ 0 < e0_. + yR0 - R > 

y Gk y z 


+ (0 z 6 Gk - Sc ( V e k )} < " yR( V B k ) " 1 


+ (R-e) {-X k R<yR(B p +B k ) + 0 Gk (-yR<; k -yR0 z -R y )-yRC k <|) k > 


+ (y R 0 y - r z > < yR (B p + B k ) -w k ? Gk > 


- A k R < yR ( B p + B k ) > 

+ [ e (0 Gk + <fr k ) + yR0 y - R z ] < yR(3 p +fy > 


+ PR (6 p + B k ) b(0 Gk + * t ) I 
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2 c, ,2 • • ' 

sinifj [ — — [ - < yR (1 + C + 9 ) - R > + eyR(R-e) ] 

K. 3. ^ K. Z X 


(R-e) 


l B k 9 * - 9 Gk'hc 9 y - 9 * < 9 Gk + +k - V - 6 Gk^k 6 k’ 


1 «k < -* ,E <6 Gk + V + 6 Gk R X ’ 


+ 0 < X. R - yR0 + R > 

x k y z 


+ 5 k < B p + V ' - BR(e Gk + V =• 


- W » e * + «V + V ‘-1* w > 


+ X. R0 - 0 < e0,„ + yR0 - R > ] 
k x x Gk y z 


+ (R-e) [ - X k R <-yR (0 Gk+ e k ) + e Gk R x > 


+ (yR0 y - R z ) < -yR (0 Gk + <|> k ) > ] ] 


2c 


+ cosij^ sim(J k [ - a 


do 


[ (R-e) yR (yR£ k + yR0 z + R y ) ] 


- (R 20 0 
3 y x 


[ -e, < 2pR (B p + B k ) - 2pK k 9 Gk > 


+ e y < -yR (0 Gk + 4> k ) > ] 


+ (R-e) [ yR0 Gk e0 y - yR(3 p +6 k ) (yR0 Gk +yR(() k ) 


+ 2 ^ R0 Gk R x (B p + B k> ] ] 


2, r (R-e ) 3 ; 2 

+ cos \p k [ 3 — 0 


i * Z f~ I 29 y < B p + V BR - ^k e Gk e y ! 


+ (R-e) [ y 2 R 2 (6 + 0 k ) 2 
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+ yR (6 

p 


. (R-e) 2 
2 

+ (R-e) 


+ S k> < (0 Gk + V + 6 Gk ( -V yR0 z )>] ] 

+ -“V i 2 

i-^Gk ( K + «kV ' "A ] 

[- U R6 Gk < e e x + (B p + B k )(-BR9 z - i y ) > 

- ~ < y 2 R 2 - 2yR R > ) ] } 

a A. 


2 

P A ab2 fi 2 { - [ ($ k - ( 6 Gk + i k )) <sin 0 Gk + 4 > k cos 0 Gk > 

- c k 0 Gk sine Gk ] 

+ (R-e) sine,,. (yR 6 - R ) 

Gk y z 

2 

- cosik [ - — R ^ (0 - 0 ) sin 0 , 

T k 2 y x Gk 

+ (R-e) [ (pR 6 k - pR( 9 Gk + 4 > k ) ) ( sin0 G k + ^ cos 0 Gk ) 


- yR£ k e Gk sin0 Gk ] ] 

- sin^ [ £=f~ (0 x + 0 y ) sin0 Gk 

• • 

- (R-e) [ (pR(6 p + 6 k ) + UR(0 Gk +<j> k )) 

(sin0 Gk + \ COS0 Gk ) 

+ yR ?k 0 Gk sin0 Gk ] ] } (5.96) 
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P A = P A a { [ 6 Gk + ^k " K + C k (6 Gk + V " C k (6 p + V 

z3k 

+ *z“ck + ^k <0 Gk + *k - + V>Gk ! 

+ [ - X k R + « (8 Gk + 8 t ) + pR9 y - E z - X k ki t 


+ e<0 Gk + *k ” V + b (0 Gk + V 1 


- (R-e) e Aj^R 


, r (R-e) ^ 

cos 4 ) k [ - 3 — 0 y 


+ [PR ( 6 p + B k ) + e Gk (-PRC k - R y - pR 0 z ) 


- “Vk + < e p + 6 k> 


+ ®Gk <- >*'k - E y - >"z> * ^k <*k - ^k) 1 
+ (R-e) [ e yR (6 p + 3 k ) + A k RyR<; k ] ] 

- slnip k [ e K 

2 

* l - uR (® G k + \) + \e Gk - pR(6 Gk+ * k )c k 

- “ R «Ak - W R<e Gk + *k - V + Vck 

- “ R< V 6 Gk + V - 6 k < 6 p +e k )+ 9 z 6 Gk > ] 

+ (R-e) [ - eyR (0 Gk + <fr k ) + A fc RyR 

- yR < e (0 Gk + 4> k ) + yR0 y - R z > 

• • 

- yR b (0 Gk + <() k ) ] ] 
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+ cos^ k sin\jj k [ yR8 


+ (R-e) [-yR<yR(£ p + + 28 Gk (-yRc k -R y -yR6 2 ) > 


- 2yR<; k 4> k 3 3 

+ cos 2 \p k [ - (R-e) y 2 R 2 ? k (£ p + 8 k )3 
2 

+ sin\ [ - yR6 x - (R-e) yR < - pR(9 Gk +\) + 2R x ® Gk > 5 > 

2 

*1 »A*» 2a2 « -• a= i i - < %- ( «Ck + +k )> 


cost P k [ (R-e) < yRB k - yR (0 Gk + > ] 

• • 

sin^ k [ (R-e) < - yR(g p +£ k ) - pR ( © Gk ' 4 '4> k > > 3 3 cos8 Gk 

(5.97) 


The aerodynamic moments at the blade root are 

• R— e 


x3k 


X 


(q A + X k C k P A " X k 8 k P A ) dX k 
x3k K X z3k K y3k K 


= p. a b n C { 


(R-e) 


[0 Gk + ^k ' B k + ^k (0 Gk + V ^k (B p + B k } 


+ 8 Ak + ^k < 0 Gk + *k _ B k > +0 Z 0 Gk ] 
+ [ - A k R + e <0 Gk + <p k ) + yR0 y - R z - A fe RC k 


+ e <9 Gk + 4-k - \) + b (6 Gk + \) 3 


+ s& f- e( -v> 


cos* t [ - 6 y 


+ &f- [ uR(6p + e k ) + e Gt (-uu k - i y - way 
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- + ? k yR A + V 


+ 0 Gk ( "^k - V yR0 z ) 


+ (^k - e k > (-yRC k ) ] 


+ (R 'f - [eyR(B p + B k ) + A k RyR^ k ] ] 


- sin* k [ ^=f- e x 


(R-e)' 


-y R (0 Gk + <f> k ) + R x 0 Gk - yR (6 Gk + <D k ) C k 


- ^ R0 2 6 Gk " yR(6 Gk + *k - 6 k> + R x 6 Gk 


y R < ^ k (e Gk + V - ^ k (e p + V + 0 z e ck > 1 


+ -— f - [ e < - y R(0 Gk + 4> k ) > + X k RyR 


- yR < e (0 Gk + <|) k ) + yR6 y - R g > - yR b (6^,+ <|> k ) ] ] 


+ cosi|J k sin^ k [ (R ~^ yR0 y 

, 2 

+ — f [-yR< yR (6 p + 3 k ) + 0 Gk (-y R ? k " R y ~ y R 6 z ) 


- yR? k 4> k > 


+ yRC k < yR (0 Gk + 4> k ) > 


- ^ R0 Gk < " yR0 z - R y >] ] 


+ cob\[ - »Vi t ce p + e t > 


+ sin 2 * k [ - (ti = | 1 nRe x 


67 



(R-e) ‘ 


[ -yR < -yR (0„,, + <pj + R x e Gk > - UR6 r1r K ] ] > 


+ i p A abVc k e k - <e Gk + + k ) > 


Gk r k 

,3 

<■ i. 

J k 


Gk x 


c°s^ k < (R ~^ [ yR ^k " yR (0 Gk + ^ ] > 


- sini|; k < (R ~ ' t - MR(B p + B k ) - yR(6 Gk + 4> k ) ] > } cos0 


Gk 


- P A abfi 2 6 k (1 + £ k 2 ) { - <l+2C k +20 z >+ 2e ] 


1 \ <(, Gk + +k - B k + V 6 Gk + V 


- c k (B p + V + e Ak =■ 


+ ^k % + V < e Gk + V ■ B k > ] 

^ # 

- ! B k < - V - 'Ak+V + >- E6 z - R z > 


• • 


+ \ R< 8 Gk + ^k ~ B k + ^k (0 Gk + V ^k (B p +B k } 


+ 0 e > 

z Gk 


• • • 

+ (-yR8 + R ) < 0_. + <j>, - 3, + C, 0-, > 

y z Gk T k k k Gk 


+ C k (B p + B k } < _X k R > + B k b (0 Gk + V 


(R-e)' 


[ " A k R < " A k R + e(0 Gk + V + yR0 y - R z > 


+ (yR0 - R ) < - A, R + e0 , + yR0 - R > 
r y z k Gk K y z 


- V> b ( 6 Ck + V 1 

cos \ t —ir 1 < i |E s. + uEO z * R y 1 
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I’ Vy ~ Wok + V 0 Gk + +k - \ > + ^ 9 Gk 6 ,] 


(R-eV 


B k < UR (S p + Sp) t e Gk <-uRt k - URG z - R y ) 


- V»5k > 


+ (X k R - uR9 y + R z > (-6 y ) ■*■ uRs k (B p + B k >' 


e Gk ‘ - ' 9 y + (6 p + V - V ’ 


+ <-* k + Bp) UR (Bp + Bp) - XpR6 y 


- h (6 ck + V >* <K + 6 k> + 6 v < e0 Gk + URB ,,-R 


y z 


+ <e z e Gk - V 6 P + V> < - «p + V " 1 
2 # 

+ — (- X k R ( UR (B p + Bp) + 0 Gk (-uRCp - URB Z - R )-URC k * k ) 


+ (uR6 y - R Z )<UR (Bp + B k ) - uR( k e Gk > 


- V yR(B P + V 


+ [ e ( e Gk + V + yR6 y - R Z ] yR (B p + V 


+ yR (g p + 3 k ) b ( 0 Gk + <D k ) ]] 


sin \ [ a 


2c do r (R-e ) 3 


[ <yR (1 + ? k+ 9 z ) - R x > + 


(R-e ) ‘ 


eyR ] 


i " * • • * ** 

■ [ae - 0., c, e - e < e„ + <*>. - e, > -e_,c,e ] 

4 k x Gk^k y x Gk k k Gk k x 


[B k < - >* R «Gk + »k> + 6 Gk E x > 


+ e x ( X k R - yR0„+Rj + <-yR(e„ 1 44>J> 


y z ' ’’k'V k' 


Gk r k' 
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+ COS^ k 


- 0_. < e0 + (0 +0, )(-pR0 - R ) + A, R0 

Gk x p k z y k x 


- 6 x (e6 ck + BR0 y - R z ) ] 


+ -a^- [- x k R < - PR <e Gt + V + e Gk \ > 


+ (yR6 y - R z )< — uR( © Gk 4> k ) >] ] 

2c ( .2 

Sin4; k [ ^ ^ ~2 [pR < pR^ k + pR0 z + R y >] 


20 0 

4 y x 


ifcfL (. 2 e x < pR(6 p + e k ) - pRc t e ok > 


+ V uR <e ck + y > ] 

2 

t [ pR 0 ok ee - p 2 rV + 6 k xe Gk+ y 


+ 2RR6 G k R * ‘ R p + V II 

+ cos 2 ip [ (R "^ e 2 
k 4 y 


!26 y (B P + y PR - pRC k e Gk e y 1 


+ [p 2 R 2 (3 + 3 k ) 2 

Z P K 


+ ^p + V < "^k (0 Gk + V + 6 Gk ( - R y^ R0 z ) >] ] 


, .2, r (R-e) „ 2 

+ sin ^ [ — 0 x 


+ iSZ f L '-^Gk ( i + ^V ' UR VJ 
2 

+ -2^- t-PR9 Gk < p9 x + <B p + 6 k )(-pR6 z - R y > > 


— < p 2 R 2 - 2pR R > ] ] } 

3 X 
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+ \ P A abV B k (1 + ^ 2 ) { - [<B k -(0 Gk +4> k )>(sin0^+4) 1 ,cose r 


Gk k Gk 


- ?k e Gk sin0 Gk ] 


+ ^=f- sin6 Gk ( ^ y - i.) 

- cost k [- (0 y - 8 x ) sine Gk 


+ [(yRB k -'UR(0 Gk +<l) k ))Csine Gk +(j) k cos0 Gk ) 


MR? k 0 G k sin 0 Gk 1 1 


- sin^ k [ ( 0 x + 0 y ) sin 0 Gk 


+ [- ( mR ( B p + 3 k ) + PR ( 0 Gk + 4 > k )) 


(sin0 Gk + ^ cos0 nv) 


Gk' 


+ pR? k 0 Gk sin 0 Gk ] ] } 


+ P A abfl 2 x A { &=f- 


• • 


[ 0 Gk + \ ~ \ + C k ( ®Gk + V ' C k (6 p + \ } 


+ 0 z 9 Gk + ^ < 8,v + - K > + 0 7®Rif ] 


Gk 


z Gk 


(R-e)' 


[- A R R + e(0 GR + 4> k ) + yR0 y - R z - A fe Re k 


+ e< 0 Gk + - B,, > + b ( 6 nh + <p^ ] 


Gk 


+ (R-e)e( -X R) 


- cos^ k [- ( -- 3 } ■ 0 y 


(R-e)' 


[ yR (B„ + B v ) + 6 Gk (-yRc v - R w - yR0,) - yR£ 


k^k 


+ c,, yR (B_ + BJ + 0 Gk (-yRc v - R, - yR 0 „) 
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+ (4> k - 3 k ) (-yR? fc ) ] 


+ (R-e) [eyR (3 p + 3 fc ) + X k RyR? k ] ] 


. , r (R-e) n 

- [ — 3— e x 


+ l-MR » Gk + V + V Gt - „r <e Gk + V c k 


- ‘‘“Ak - >® < % + ♦k - 6 k > + R Ak 


~ VR< h (0 Gk + *k> - S k «p + V + 0 z 0 Gk > 1 


+ (R-e) [ e < - yR (6^ + ^ > + X fc RyR 


- yR < e (0 Gk + 4> k ) + yR6 y - R z > 


" yR b (0 Gk + V 1 3 


+ cos^ k sini^ k [ ^ R ^ — yR0 y 


+ (R-e) [-yR < yR(3 p + B fc ) + 6 Gk (-yR? k -R y -yR0 z ) 


- t«kV 


+ ^k < ^ R (0 Gk + +k ) " 


- yR0 Gk < -yR0 z - R y > ] ] 


+ cos\ [- (R-e) y 2 R 2 C k (p + B fc ) ] 


+ sin 2 ^ k [- yR6 x 


+ (R-e) [-yR < - yR(0 Gk + <^ fc ) + > - yR0 Gk R x 1 1 > 
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. 1 , 2^2 , b s f (R-e) / * , * % 

+ 2 P A a b a (X A - 2> { 2 <3 k ~ (0 Gk + V > 


- cosi|J k < (R-e) [pR3 k - URO Gk + 4> k ) ] > 


• • 


- sini|> k < (R-e) [-yR(B p +B k )-yR(0 Gk +(j) k )]>} 


- V «Gk + V « 2 1 V ^ < R - e > 




(R-e) ‘ 


2 ' Gk Y k' 


(R-e) e 


- (R-e) < yR? k cos^ k + yR sin^ k > } 

(5.98) 


y3k 


/•R-e 

^ r, 


(c *a ^ + X A p A ' (x k + V p a > dx k 

y3k x3k z3k 


Neglecting higher order terms 


y3k 


-x 


R-e 


- x k P A dx k 
k z3k K 


• • 


p A abn 2 ( [0 Gk + - B k '+ 


+ 0 z 0 Gk + ^k^Gk+W + 6 z e Gk ] 


+ i5=|2_ [- x k R + e(0 Gk + + k ) + yR0 y - R z - x k R c k 


+ e(0 Gk + *k " V + b(0 Gk + V 1 


+ e (-x k R) 


(R-e) 


4 . 


- cos^ [- 0 
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(R-e) 


+ 3 ll,R(|i p + V + 6 Gk ''^k - R y - >*V 


' ^k \ + ^k RR lB p + V 


+ e Gk (- U Rt k - R y - pR0 z ) 


- ‘"Sc »k - B k» 


+ ■ ( * r f ) I euR (B + 6 k ) + » k «kRC k 1 1 


- sini)J I BR ? ^ — S 
k 4 x 


(R-e) 


+ l-wR(e Gk + » k ) + v Gk - pr (e ck + * k ) t k 


-pR8 z e Gk - uR(e ck + * k - B k > + y ck 


- BE ‘ H <6 Gk + *k> - V B p + B k> * 1 

2 

+ te < - pR (6 Gk + 4> k ) > + X k RpR 


pR < e (0„. +<}>.) + pR0 - R > 
Gk k y z 


- PR b (0 Gk + 4» k ) ]] 


+ cos^ k sin^ k [ — R pR6y 

+ I-yR < PR(B p + 3 k > + 8 Gk (-pRC k - R y - PR6 z ) 


- ^kV 


+ ^k^ R (9 Gk + *k> “ ^Gk " - **». ~ R y > 1] 


+ cos 2 ty k [ - pV^ (B p + B k > 1 


74 



+ sin 2 i}j k [ - (R ~3 ' UR6 x 


+ I' " R •= ' ^ < B Gk + V + V>Gk =■ - '“cA 111 

< K - «gr + v> 

- cos^ k < (R ~^ ' [ yR0 k - yR (9 Gk + 4> k ) ] > 


• • 


sln V [ - pR (e p + 3 k } " pR (e Gk + K )]>} COS0 Gk 


z3k 


/.R-e 

Jo 


(5.99) 


(q 4 + (x k + V P A " W P A ) dX k 

z3k y3k x3k K 


Neglecting higher order terms 


.R-e 


z3k 


/ (x k P A ' x k ? k P A „ > dx k 

j o 


y3k 


x3k 


where 


x3k 


' 5 k P A - 6 k P A = "Sc P A \ P A 


y5k 


z5k 


y3k 


z3k 


thus 


z3k 


= r e ^ v + ^ ^ ^ 

^ 3 

= p . a bfi 2 { - — [ (R ~^ < 1 + 2C + 20 > + (R ~! -— 2e] 

A a . -9 k z 3 


• • 


^V- < B k ‘ 0 Gk + *k - 6 k + C k «Gk + V ‘WV 


+ e e„, > 

z Gk 


+ h < b p + V <e Gk + h - B k >] 


i5= | L - IB k < - V + » (9 Gk + *k> + - E z > 
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+ V * 0 Gk + +k " 6 k + S <e Gk + V ' \ (E p +e k > + e z ®Gk > 


+ <-«*, + V •= 6 Gk + *k - 6 k + VGk > 


+ S k<*p + 6 k )< - X k R> +6 k b (e Gk + V ] 

2 

+ (- k k k < - V + el6 Gk + V + *"y - V 


+ (yiR0 - R ) < -X, R + e0-, + yR0 - R > 
y z k Gk y z 


- V b<6 Gk + V 1 

2c , >. v 3 • 

- cos4/ k [ [ 1 ~y < yRC k + UR0 z + R y > 1 


(R-e) 


- [-M.. - 5,1.^,. + e„ < 0 nu + - B, > + L6 P1 A, ] 


4 L V y k x Gk T y Gk ^k \Gky 

i\ •= ee < B P + v + v-^k - - y-w 


+ <V - bE9 y + E k )< - e y ) + S k (B p + V EE < B p + V 


-0 Gk < -a0 j+ <6 p + P k )(liR - R x )> 


+ c-* k + B k ) PR (B p + B k ) - >' k BO y 


- 'k (6 Gk + V EE < 6 p + V 


+ 0 (e0 , + yR0 - R ) 

y Gk y z 


+ (9 2 9 Gk - 9 k « P + y> < - ue(b p + v* 

2 • 

+ I- X k R < PR (B p + B k ) + 8 Gk (-pRC k - PRB Z - » y ) 


- ‘"SV 
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+ <>* R6 y - R Z> < •* «p + \> - •"« k e Gk > 


X k KpR ( 6 p + 6 k ) 


+ (« <e Gk + * k ) + m y - R z ) pr <s p + B k ) 


+ pR <B p + 6 k ) b (9 Gk + \) 1 1 

- sl ”* k t t ss= f^ < 1' R ( 1 + + V - 

,,4«. *• • i • 

- - -J 1 - [g, 0 - e„, c e - 0 < e„, + <)). - 8, > - e„, t. 0 

4 k x Gk k y x Gk T k k Gk k 2 


^T" ' 9 k ‘ - RR ">G k + V + 9 Gk V 


* \ < V - l* Re y + R z ' + h «p + 8 k> < - RR(6 Gk + V > 


9 0 k ' ' 6 X + “p + V^z - V " + X k R9 * 


- 6 k <et, Gk + RRe y - R z> 1 


+ - & f“ 1 - X k R < - RR «Ck + V + W 


+ (yR0 - R ) < - yR (0„, + 4>, ) > ] ; 

y z IjK K 

2c , s2 

+ cosij;, siniK [- — — ^ , [ yR < yR£, + yR0 + R > ] 

k. k. el t- K z y 

- fR ~*' J 29 f 
9 y 


-&f- [ - 29 x < PR (B p + e k ) - RRC k 6 Gk > 


+ 9 y * -" R (9 Gk + V 0 1 


eyR ] 

] 
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+ [ MR6 Gk e9 y - pV CB p + 6 t ) « Gk + * k ) 


- 2 “ R6 Gk R x <B p + V 1 1 

+ =“\ <■ K 2 - ■ a f- t2 V e p +6 k )UE ' '“VW, 1 

+ < fc f>- 2 [1J 2 R 2 (6 + B > 2 

L p K. 

+ + V •= - «k (e Gk - V + 6 Gk <- \ - "*,» > 1 1 

- »*»\ [ K 2 + [ -' JR6 Gk ‘W^W 

( o_ e )2 • * 

+ [ -UR0 p , < ee + (6 + 6, ) (-UR0 - R ) > 

2 Gk x p k z y 


“do . .22 


< y R - 2yR R > ] ] } 


3 

I P A a b ^ 2 { - l< h’ (0 Gk + V > (sin0 Gk + \ COS0 Gk ) 


C k 9 Gk Sin9 Gk ] 


+ sin0 Gk ( P R0 y - V 

- C0S *k [ - SS lh (0 y - V sin0 Gk 


+ [(yR0 k - yR (6 Gk + ^))(sin0 Gk + <J> k cos0 Gk ) 


- yR^ k 0 Gk sin0 Gk ] 

- si# k [ (0 x + 0 y ) sin0 Gk 

2 

+ ^=f*- [ - (PR (B p + B k ) + PR W Gk + i k )) 

(sin6 Gk + h COS0 Gk ) 
+ PR? k 0 Gk Sin0 Gk ] J } 
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+ p A .i»fl 2 e k c k ( <% + * k - V + <-\w 

- cosip k [ (R ~^ yR (B + B k )] 

3 • 2 

- sintj^ [ - ( - R ~ 3 } [- 2yR (0 Gk + «j> k ) + yRB fc ] + (R ~| ) -- ■ X fc RyR ] 

2 

+ cos^ k sii# k [ - Y~ ^ 2r2 (3 + B k ) ] 

+ sl »\ 1 15= f^ • j2r2 «Gk + V 1 > 


(5.100) 


The orders of magnitude of the leading terms in the various aerodynamic loads 
are given for convenience 



= 

0(e 3/2 

x3k 


= 

0(e 3/2 

y3k 
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o( E 1/2 

z3k 


= 

0(e 3/2 
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= 

0(c 1/2 

y3k 


= 

0(e 3/2 

z3k 



) 

) 

) 

) 

) 
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Damping Force 

Instead of assuming a distributed damping force representing the structural 
damping of the blade, one can assume a damping force proportional to the 
velocity and it is provided at the root of the blade. It can be written as 
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Flap 

V ' 

+ e k 


(5.101) 



b F 


Lead-lag 


- n k 

b L 

(5.102) 

Torsion 

\3k 

- a \ 

8 s x 

(5.103) 


These damping forces are proportional to the rate of change in angles in flap, 
lead-lag and torsion. Reason for this type of assumption is that in the 
present study the blade is modeled as a rigid blade with root springs. 

Equations of Motion 

Using the moment equilibrium condition at the blade root, the equations of 
motion for the k*-* 1 blade can be written as follows. 


Flap: The flap equation is 

M + Q + Q + Q =0 
k y3k y3k y3k 


(5.104) 


Lead-lag: The lead-lag equation is 


Mr + Qi + Q A + Q D “ 0 

^k z3k z3k z3k 


(5.105) 


Torsion: The torsion equation is 


M * + Qi + Q a + q d 
9 k x3k x3k x3k 


(5.106) 


The elastic restoring moments due to the root springs are given in Appendix. 
The results are summarized below. 




<MkV 




A 


{ K g + Ri (K 5 - K p ) sin 2 0 Gk ) 
a k +4 ’k B k ) 

+ R x (K e - Kg) sin0 Gk cos0 Gk 



R ! ( \ " K B> sin20 Gk } 

<VVk } 

-T R x (K^ - K g ) sin0 Gk cos©^ 
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\ " ‘ K +1 ^ ^ ’ " V k 


where 


k r k r 
b h 3 b 

k r +k r 


K K 
K +K 


1/K. fi - 1/K 
g B ? B 


(1/K fl + 1/K„ ) - (1/K + 1/K ) 

g B 3 H S C H 


A' = 1 + R l (1-f^) sin 2 6 Gk (K^ - K g ) 2 /K ? Kg 


For an articulated blade 

K e - K c -v - 0 and V 


K , k a 
^B 
K a +K , 
*0 


For a hingeless blade 


K, =0 and K 


K . 
*0 


'+1 K <j> +K «f> 
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6. EQUATIONS OF MOTION OF THE SUPPORTING STRUCTURE 


6.1 General 

The supporting structure, to which the envelope and the rotors of the HHLA 
vehicle modeled in this study (Fig. 2) are attached, is assumed to be flexible. 
The degrees of freedom associated with the supporting structure consist there- 
fore of both rigid body and flexible degrees of freedom. The structure is 
idealized as a free-free beam and is represented by relatively few number of 
bending and torsional free vibration modes. The derivation of the equations of 
motion consists of two parts. 

1) One part representing the rigid body degrees of freedom. 

2) A second part representing the flexible modes of the structure. 

In dealing with these two contributions to the total motion it is assumed 
that rigid body motions occur first. After establishing a perturbed state, con- 
sisting of rigid body translation and rotation, the elastic deformations are 
superposed. Thus the equations representing the rigid body degrees of freedom 
are written in S-system and the equations of motion representing the elastic 
structure are written in Sl-system, see Figs. 8 and 9. To write the equations, 
various external loads have to be identified and properly transformed to the 
corresponding coordinate systems. 

6.2 Loads 

The various external loads acting on the model are illustrated in Fig. 12, 

1. The rotor loads at hub center 0 U1 and 0,, K . 

HI H2 

2 . Aerodynamic loads on the envelope acting at the center of buoyancy Op, 
and a static buoyant lift. 

3. Thruster loads acting at and 

4. Gravity loads acting at the respective center of mass of underslung 
weight, envelope weight, fuselages weight and the supporting structure 
weight. The respective center of mass locations are denoted by Ojjpj, 

°EN’ °F1> °F2’ °S* 

6.2.1 Rotor Loads 

The loads acting on the rotor blades are inertia and the aerodynamic loads. 
These distributed loads are transformed into forces and moments acting at the 
blade root. The moments are balanced by root springs and dampers. The loads 
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on the k*^ blade are derived in Chapter 5 and are given in 3k coordinate system. 
They are given in equations (5 . 89) — (5 .100) . Net combined loads acting at the 
root of the blade are 


P = p + p 

r 3k I3k A3k 


Q 3k Q l3k + Q A3k + Q D3k 


( 6 . 1 ) 

( 6 . 2 ) 


These loads are transformed into loads acting parallel to 2k system, using the 
following transformation. 


' 2k3 


[t 23 ] 


‘3k 


(6.3) 


^2k3 ^ T 23^ Q 3k 


(6.4) 


where P^^ anc * ^2k3 are f° rces an ^ moments acting at the blade root (i.e. 
origin of the 3k system) whose components are given in 2k system. Transforming 
these loads into loads acting at the origin of the 2k system (i.e. hub center, 
either 0 Rl or 0^) 


2k = 

P 2k3 





(6.5) 

'2k = 

^2k3 + 

(e 

e x2k X P 2k3^ 



(6.6) 

2k = 

(P x3k " 


i P z3k )e x2k + 

A 

P y3k®y2k + 

(P z 

3k +P p P x3k )e z2k 







(6.7) 

2k = 

(Q x3k “ 

3 

P 

i Q z3k )e x2k + 

Q y3k e y2k + 

( ^z 

3k +P p Q x3k )e z2k 


+ 

[- 

■e (P ,, + 3 
z 3k p 

P x3k ) e y2k 

+ e 

P e ] (6.8) 

y3k z2k 


Again transforming these loads to 1 system 


P lk [T 12 1 P 2k 


(6.9) 


^lk [T 12- 1 Q 2k 


( 6 . 10 ) 
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Expanding 


P 


lk 


Q 


lk 


= [cos^ k < P x3k - g p 

P z3k > 

sirajj 

k P y3k 

1 e xl 




+ [s 

ln \ 

< P x3k 

6 p P z3k 

> + cos^ k Py3k ] a yl 



+ [P 

oi + $ I 

z3k p 

‘x3k^ 

A 

e zl 


(6.11) 

1— 1 
X 

< 0) 

II 

[ COSl|l k 

< Q x3k 

- 6 

P 

Q z3k 

> - 

sint(i k < 

Q y3k 

- e(P z3k +P p P x3k )>] 

A 

+ e . 

yi 

[ sin^ k 

< Q x3k 

- 0 

P 

Q z3k 

> + 

cosilt < 
k 

Q y3k' 

' e(P z3k + V’x3k )>] 


+ S zl [ «z3k + 6 p «x3k + e P y3k ] 


These loads are acting origin of the 1 system i.e. either at 0^ or 0^ 2 de- 
pending on the rotor system. Total load due to all the blades in the two 
rotor systems are 


N 

Rotor System 1 ^1H1 = S ^lk (6.13) 

k=l 


Q 1H1 £ ^lk 
k=l 

N 

Rotq^SjrstegLZ P m = £ P lk 


(6.14) 


(6.15) 


N 

Q 1H2 = ^ Q lk 

k=l 


(6.16) 


The summation is over the number of blades in the individual rotor systems. As 
it was pointed out earlier that P^ and Q^ k will be different for different 
blades and also for different rotor systems. 

Transforming the loads acting at 0-^ to the point 0p 3 (C.G. of fuselage Fl) 
and that at 0g 2 to the point 0p 2 (C.G. of fuselage F2) , 
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(6.17) 


P = P 

r lFl 1H1 


Q 1F1 Q 1H1 + (h 2 e zsl x P 1H1 5 


Q 1U , + P e . + h„ P e ] 

1H1 2 y 1H1 xsl 2 x m ysl J 


(6.18) 


where P is the y-component of the vector P. and P is the x-component 

ym x ihi 

of the vector It should be noted that according to our initial assumption 

the 1 system and the SI system are assumed to be parallel. Similarly for the 
second rotor system 

P lF2 = P 1H2 (6.19) 


Q 1F2 Q 1H2 + (h 2 e zsl X P 1H2^ 


= Q lu0 + [-h P e . + h P e ] 
1H2 2 y 1H2 xsl 2 x 1H2 ysl 


( 6 . 20 ) 


The loads presented in equations (6.17)-(6.20) are acting at the points 0^ and 
0 . The components of these loads are given in 1 system (SI system) , both 

the 1 system and SI system are parallel systems. The components of these loads 
will be used in writing the equations of motion representing the elastic modes 
of the supporting structure. Writing these loads in terms of loads in the 3k 
system 


N 

P 1F1 = £ xl fcos \ < P x 3k " P p P z3k > " Sln \ P y3k^ 


+ e yl [sin^ < P x3k - 8 p P z3k > + cos^ k P y3k ] 


+ e zl [P z3k + » p ^k 1 } 


( 6 . 21 ) 


'1F1 


■E ( e [cosij; < Q„„, - > - sini^ < Q, 

k=l 


xl l “““ r k ' ^x3k M p x z3k ' ' “*“^k ' x y3k 


e (P z3k + g p P x3k ):> 


-h 2 (sin^ k < P x3k - 8^ P z3k > + co S1 p k P y3k )] 
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+ e yl [sin* k < Q x3k - Q z3k > + cos^ < Q y3k 


e(P z3k + g p P x3k^ > 


+ h 2 (cosiJ; k < P x3k - 3 P z3k : 


Sin \ P y3k ) ] 


+ e 


zl 


"W + »p ‘W + e P y3k 1 > 


( 6 . 22 ) 


The expressions P^^ and represent the forces and moments due to rotor 

systems Rl. These loads act at the point 0^ on the supporting structure. 
Similarly, the loads due to the rotor system R2 are P^ 2 and ^1 f 2‘ These loads 
act at this point 0„. on the supporting structure. Since the rotor loads are 
derived for a general rotor system, the expressions for these loads given in 
equations (6.21) and (6.22) are valid for both rotor systems Rl and R2 . How- 
ever the components P x3k> P y3k> P z3k> Q x3k , Q y3k , Q z3k could be different for 
the two rotor systems depending on the operating conditions and various other 
rotor parameters. 

The rotor leads ^1F1’ P 1F2’ ^1F2 t> e usecl i- n writing the equations 

of motion of the vehicle and the supporting structure. Hence, it is convenient 

to refer these loads in the body fixed Si coordinate system. From equations 

(6.21) and (6.22) it can be seen that these rotor loads are given as components 

along the hub fixed 1 system. By definition, the hub fixed 1 system and the 

body fixed SI system are parallel. Therefore in subsequent parts of this report 

these rotor loads will be referred as P c for P nT . 1( Q„ for Q, P„ for 
_ b lFl iF1 S 1F1 S 1F2 

P, „ and Q„ for Q 1lr0 . Again these loads are transformed into components along 
Xr Z xr z 

the S system acting at point Og on the structure. Finally the individual rotor 
loads are added together to get the total rotor loads. 

The total rotor loads at Og are (components along S system) 


P SR ^ T S,S1-* P S1F1 + t T S,Sl ] P S1F2 


(6.23) 
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Q SR * T S,S1^ ^ Q S1F1 ~ ^Fl e xSl x ^SlFl^ 

+ ^ T S,S1^ ^SlF2 + £ F2 6 xS1 X ^S1F2 ^ ( ~ 6 - U '> 

where [Tg g^] is the transformation matrix relating the body fixed noninertial 
SI system to inertial S system. 


[T 


S,S1 J 


0 0-0 
y x z 


0 0+0 

z x y 


0 0-0 
z y x 


L -0 


where 0^, 0^, 0^ represent the Euler angles for rigid body rotation, 

6.2.2 Aerodynamic Loads Due to the Envelope 

The aerodynamic loads of the envelope are the buoyancy loads acting along 
the body axes. These are denoted respectively by P^ and Q EN * 

P„„ =P e ,+P e.+P e ., + P S e 
EN x xsl y ysl z zsl z zs 

= (P - 0 P S ) e 1 + (P +0P S )e . + (P +P S )e 
x y z xsl y x z ysl z z zs 

(6.25) 


Qt-n' - Q ^ . + Q e - + Q e 

Els x xsl y ysl z zsl 


(6.26) 


These loads act at the center of buoyancy 0^. The components of the dynamic 
loads are defined as [Ref. 7] 

Forces P = C q V 2/3 , P = C q V 2/3 , P " C q V 2/3 (6.27) 

xx yy zz 


Moments 


Q = C. q V , Q = C q V ,Q=CqV (6.28) 

n m z n 


'‘x £ ‘ ' 'y m 

Where V is the volume of the envelope 
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q is the dynamic pressure 


and 


P is the static lift on the envelope acting along Z axis 
z s 

c , c„ are coefficients 

x, y, z x., m, n 


Transforming these loads to the point 0 on the supporting structure 


P S1EN P EN 


(6.29) 


Q S1EN ^EN + (h 3 e zsl X P EN ) 


Expanding 


Q EN + (P - + e -. p !) e__, + h, (P_ - e.P!)e._ v ] (6.30) 


-> y x z xsl 3 x y z ysl 

2/3 


• Q1ffw = (C qV 2/3 -0P S )e + (C q V^ ,J + 0 P " ) e , 

S1EN x n y z xsl y M x z ysl 


+ (C q V*' J + P s ) e n (6.31) 
z z zsl 


<SlEN 


(C q V - h. C q V 2/3 - h,0 P S ) e . 
H 3 y 3 x z xsl 


+ (C q V + h, C q V 2/3 - 0 P S )e . 
m 3 x ^ y z ysl 

+ C q V e . (6.32) 

n zsl 


Writing these loads. 

along the 

S system, acting at 0g 


p 

SEN 

^ T S,S1^ 

P S1EN 

(6.33) 

Q SEN 

^ T S,S1^ 

Q S1EN 

(6.34) 


6.2.3 Thruster Loads 

Two thrusters simulated by thruster forces are assumed to act at the C.G. 
of the fuselages Fl and F2 i.e. 0^ and 0^- They are 


A 

P = - p e ' 
TFl T xsl 


and 


P = - P e 
TF2 T xsl 


(6.35) 


Transforming these forces to the C.G. of the structure 0 and combining them 


SIT 


2 P m e . 
T xsl 


(6.36) 
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in the S coordinate system 


P ST ^ T S,S1^ 


SIT 


(6.37) 


6.2.4 Gravity Loads 

The gravity loads are due to the various masses which constitute the 

A 

HHLA model. There loads act along ~ e zs direction, at the respective centers 
of mass. The various contributions are: 

Fuselage 

The gravity loads on the two fuselages act at 0_, and 0 „. These are 

r X F z 


p = - W e 
WFl Fl zs 


(6.38) 


" WF2 


W r , e 
F2 zs 


(6.39) 


Transforming these loads to the point Og on the structure and adding them 


P SWF - (W F1 + W F2 ) e zs 


(6.40) 


Q SWF ~ ^Fl e xsl X P WF1 + £ F2 e xsl X P WF2 

where P gWF is the gravity force 

Q SWF is the moment at Og due to the gravity forces. 

Qcinr = (“X>i e „ ~ 6 e +0 A,,, e ) x (-W e ) 

X SWF Fl xs z Fl ys y Fl zs Fl zs 

+ (£ „ e +0 £ e - 0 £ „ e ) x (-W „ e ) 
F2 xs z F2 ys y F2 zs F2 zs 


<Vf1 W F1 - "V 


+ <-t n H pl + » F2 ) e yB 


(6.41) 


Underslung Weight 

It is assumed that the underslung mass is rigidly attached to the support- 
ing structure. The gravity force on this mass acts at its C.G.Oy^ at a 

A 

distance -h.e . from 0„. 

1 zsl S 
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The loads due to this mass at Og (in S system) are 


P = - W e 
SWUN UN zs 


(6.42) 


Q SWUN h l e zsl X ^ W UN e zs ) 


= (- h (0 0 + 0 ) e - h 1 (0 0 - 0 )e -he )x(-W e ) 

1 z x y xs 1 z y x ys 1 zs UN zs 


= h. (0 6 - 0 ) W e - h. (0 0 + 6 ) W e (6.43) 

1 z y x UN xs 1 z x y UN ys 


Envelope Weight 

The gravity load on the envelope acting at its C.G., 0_», is -W_,, e . It 
° EN EN zs 

is located at a distance h^ e ^ from Og. Transforming this load to the 
point Og on the structure. 


P = - W e 

SWEN EN zs 


(6 .44) 


= h. e x(-W_ x , e ) 
SWEN 4 zsl EN zs 


[h. (0 0 + 0)e + h. (0 0 - 0)e + h. e ]x 

4 z x y xs 4 z y x ys 4 zs 


( '"en S s > 


■-»,(•! - 0 ) e n, (1 1 +6 )M e (6.45) 

4 z y x EN xs 4 z x y EN ys 

Supporting Structure Weight 

The gravity force on the supporting structure is acting at its C.G., Og, 

in -e direction, it is given by 
zs 


P = -W e 
SWS S zs 


(6.46) 


Q SWS h 5 6 xsl X * W s £ zs^ 


(h_ e +0he - 0 h. e ) x (-W e ) 
5 xs z 5 ys y 5 zs S zs 


(-0 h W + h, W e ) 
z 5 S xs 5 S ys 


(6.47) 
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6.3 Rigid Body Equations of Motion 

Using the various loads derived in S coordinate system, the rigid body 
translational and rotational equations of motion can be written in the S 
system. 

Let the rigid body perturbational translational motion of the point 0 g 
be 

R = R e + R e + R e (6.48) 

s xs xs ys ys zs zs 

Then the translational equations of motion become 


g «s P SR + P SEN + P ST + P SWF + P SWUN + P SWEN + P SWS (6-^9) 

where W = W £N + Wy N + W pl + W p2 + W g 
Wg is the supporting structure weight. 


The equations of motion for the rotational degrees of freedom are 


at- [I Qqd Qci.n? + 


+ Qc 


+ 


(6.50) 


where [1] = [Tg^ g ] [ X g + I Fl + I p2 + I £N + I UN 1 t T Sl>S ^ 

and o)= (0 - 0 0 ) e + (0 + 0 0 ) e + (0 - 0 0 ) e 1 

x zyxs y zxys z yx zs J 

The individual inertia tensors are given by 
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6. A Equations of Motion for 

Elastic 

Modes 






The flexible supporting structure is assumed to be a free-free beam 
idealized by a small number of modes of vibration. The structure can have 
bending vibration in the two planes and a torsional vibration about its 
longitudinal axis. These elastic modes are assumed to occur in SI coor- 
dinate system, and therefore the equations of motion are written in SI system. 
Furthermore the structure is assumed to vibrate about its equilibrium position. 
Bending in x gl - y ^ plane [Fig. 9] 

The elastic deformation at any point on the structure due to the bending 
deformation is written as 


v 

s 




NM 


1=1 


n i U sl ) 5 i (t) 


(6.51) 
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where is the mode shape, is the generalized coordinate or normal coordinate 
and NM is the number of modes. 

The equations of motion written, in normal coordinate, are 

['MsK?} + fKj {£} = {Q} (6.52) 

Where { Q } represents the generalized forces 

f'Mv] is the generalized mass 
f'Ks) is the generalized stiffness 

It is assumed that these vibration modes are small compared to the rigid body 
perturbational motion. So there is no coupling between the vibration modes. The 
generalized forces are (force at a location x modal displacement at the same 
location or moment x slope of the modal displacement curve) 

Q i = P ySlFl n i (il Fl ) + Q zS 1F1 n i,x (5 'Fl ) + P ySlWFl n i (£ Fl ) 

+ P ySlF2 n i (A F2 ) + Q zS 1F2 n i,x (£ F2 ) + P ySlWF2 n i (K 'F2 ) 

+ P ySlEN n i ( °S ) + Q zS 1EN n i,x ( °S ) + [P ySlWEN + P ySlWUN ] T1 i ( °S ) 

+J F p ySlws n ± (x) dx 
" £ fi 

Terms with suffix 'S1W' refer to gravity loads referred in SI system. A typical 
element of the mass matrix and stiffness matrix are 


m f £p2 / ^ , , W F1 2 

M ii = J n i m s (x) V X + T" 


l <hl> + "l,* «F1> 


+ V r, 2 «„> + I n i x «F2> 

g 


+ -f n i ( V + X ZZ EN p i, 2 x <V 


* ~f ^ ( V + T ZZ un n i,x <V 


2 f ^F2 

and K. . = oj. / n. m (x) n . dx 

n i J 'is 'i 


93 



where is the mode natural frequency of the supporting structure 
m g is the mass of the structure/unit length 
Bending x g ^ - z plane 

The elastic deformation at any point is given by 

NW1 

w s (x si ) = S n.(x ) £ (t) 
s s i=NM+l 1 sx 1 


(6.53) 


The form equation of motion is similar to equation (6.46). The generalized 
forces are 

Q i = P zSlFl n i (S 'Fl ) + Q ySlFl n i,x (5 'F1 ) + Q ySlWUN n i,x (0 S ) 

+ P zSlF2 'WP + Q ySlF2 n i,x (S 'F2 ) + Q ySlWEN n i,x (C V 


+ P zSlEN n i (0 S ) + Q ySlEN n i,x (( V 
Typical elements of the mass and stiffness matrix are 


M ii = i : 1 n i m s (x)ri i 

F2 


W 


dx + 


W 


FI 

g 

! F2 


g 


W 


W. 


UN 


and 


2 f £ Fl 

C. . = 0). / U . 

11 1 0 1 
F2 


m (x) u. dx 

S 1 


n i ( V + I yyFl n i,x (£ Fl ) 
"i ( *F2> + I yyF2 *i .VV 


f "i <°S> + P yyEN "i.^V 

W 

, EN , 2 2 ... 

+ — ‘ "i,* <°S> 

"i <°S> + 'yym "l.VV 

w 

+ T“ h l n i,x ( V 


94 



Torsion 


The torsional vibration of the supporting structure is about x , axis 

si 

[Fig. 9]. The twist at any section due to the torsional deformation is given 
by the normal modes n^( x ) as 


NM2 

*8 (X S1 > ’ . £ 5 i <0 

x=NMl+l 


(6.54) 


The torsional vibration equation in normal mode is 

• • 

(U + TKJ U) = {q} 


(6.55) 


The generalized force Q^'s are 


Q i Q xS 1F1 n i (Jl Fl ) 


+ Q xS 1F2 n i^F2’ ) + Q xS 1EN n i 
+ ^xSlWEN n i + Q xslWTJN n i^°S^ 


Typical members of [Tj and f'Kj matrices are 


I. . 

11 


C 2 ", 

*F1 


and 


K. . 

11 


= 0 ). 


i s ( x) n i dx + i xxF1 ^ U F1 > + I xxF2 ^ (^ F2 ) 


+ T XXEN n i (0 S } + r xxUN ^i (0 S ) 

w w 

+ ~T h 4 n i + T" 1 n i (0 S } 


/ «, •_ 
F2 

FI 


I (x) ri . dx 
s x 


where I is the moment of inertia of the structure per unit length about its 
longitudinal axis . 
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7 . CONCLUDING REMARKS 


A complete set of dynamical equations of motion for a simple model of 
HHLA were derived in this report. These equations can be used to study the 
stability of HHLA and to obtain the various response quantities at different 
stations on the vehicle. For convenience, the equation numbers are summarized 
below with the physical degree of freedom which it represents. 


Blade equations 


Flap Equation (5.104) 
Lead-Lag Equation (5.105) 
Torsion Equation (5.106) 


sorting Structure equations 


Rigid body translation Equation (6.49) 

Rigid body rotation Equation (6.50) 

Bending in x s ^ - y s ^ plane Equation (6.52) 

Bending in x ^-z plane Equation (6.52) 

Torsion about axis Equation (6.55) 

The coupled rotor/body equations of motion which have been derived in this 
study have considerable versatility and can be used to model a number of diverse 
rotary-wing configurations, which are listed below: 

(a) Isolated rotor blade aeroelastic stability. 

(b) Coupled rotor/ fuselage dynamics for a single rotor. 

(c) Response cyclic, collective and higher harmonic control inputs. 

(d) Stability analysis of a tandem rotor system connected by a flexible 
structure . 

(e) Dynamics, aeroelasticity, and aeroelastic response of a Hybrid Heavy 
Lift Airship. 

Depending on the type of system which one intends to analyze the complete 
equations presented in this document have to be simplified to fit the specific 
application. 

In a sequel to this report entitled "Aeroelastic Effects in Multirotor 
Vehicles, Part II: Method of Solution and Results Illustrating Coupled Rotor/ 

Body Aeromechanical Stability", two separate coupled rotor/body problems are 
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solved with considerable detail. In the first case the equations are used 
to predict the aeromechanical stability problem of a single rotor helicopter 
in ground resonance, including the effect of the aerodynamic forces. For 
this case high quality experimental results are available, and the agreement 
between theory and test was found to be quite good. In the second case, the 
stability of a simplified model vehicle (Fig. 2) representing an HHLA type 
vehicle in hover is analyzed, and the basic aeroelastic characteristics of 
such a vehicle are obtained. 

The various details of the solution such as: evaluation of the equilibrium 

position, stability equations in multiblade coordinates and appropriate methods 
of solution are given in the second report, which constitutes a sequel to the 
present report. 
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Figure 2. Sketch Showing the Main Ingredients of the Aeroelastic Model 



102 



Figure 3. 


Blade Cross-Section Configuration 
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Figure 6. 
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Figure 8. Inertial S-Coordinate System (Coincident 

with Body Fixed Si-System in Unperturbed State) 


Figure 9. Body Fixed Si-Coordinate System (.Obtained After Rigid Body Perturbationai 

Motion) Elastic Deformations of the Supporting Structure Occur in this Si-System 



y 



Figure 10. Geometry for Oscillating Airfoil in Pulsating Flow 
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Figure 11. Relative Flow Velocities 

Lq - Circulatory Lift Normal to the Resultant Flow 
L^jq - Noncirculatory Lift Normal to the Blade Chord 
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Appendix A: Equivalent Structural Model for 

Offset Hinged Spring Restrained Blades 

A. 1 General 

In this Appendix the expressions for the elastic restoring moments acting at 
the hinge of a spring restrained blade are obtained. The root springs are re- 
presentative of blade structural flexibility or they can represent flextures built 
into the blade, thus the model simulates the elastic properties of the configuration 
shown in Fig. A.l. The various expressions for the elastic restoring moments in 
flap, lead-lag and torsion, respectively, are derived. Subsequently these expres- 
sions are compared to similar equations obtained in previous studies. In this 
comparison both elastic hingeless blades as well as spring restrained equivalent 
models are considered. Furthermore it should be noted that the main advantage in 
using this simple model for an HHLA type vehicle consists of the capability of the 
model, to capture the essential behavior of both hingeless and articulated rotor 
configurations. 

Ormiston and Hodges [Ref. 13] have derived one of the first models of this type, 
however Ref . 13 was restricted to the equivalent model of a hingeless spring re- 
strained blade having only flap and lag degrees of freedom. The present model 
represents an extension of Ref. 13, to the case where both torsional blade flexi- 
bility as well as pitch link flexibility are incorporated in the blade model. 

Peters [Ref. 14] has also derived a flap-lag model similar to Ormiston and 
Hodges. In both studies, Ref s . 13 and 14, the flexibility of the hub has been also 
considered. A careful study of these two References reveals some descrepancies , 
which are associated primarily with the hub flexibility . The expressions for 
hub flexibility are also derived for our model since it was felt that these could be 
useful in some potential application. The final moment expressions obtained in this 
study are capable of representing the coupled flap-lag-torsional motion of both 
hingeless and articulated blades. 

This Appendix is divided -into three sections. In the first section the moment 
expressions, excluding hub flexibility, are derived for a spring restrained blade. 

The second section treats the special form of these equations for the case of an 
articulated blade. Finally the last section presents the moment equations for 
the blade including hub flexibility. 
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A. 2 Moment Equations for a Hingeless Blade 

The flexible elastic blade is modeled as a rigid blade with three orthogonal 
springs located at the root of the blade. These springs represent the flexibility 
of the blade in flap, lead-lag and torsion, respectively. In addition to these 
three springs, a torsional spring, in board of these three springs, is introduced 
to represent the control system or pitch link flexibility, see Fig. Al. It is 
assumed that the orientation of these springs does not change as the blade under- 
goes reorientation due to flap, lag and torsional motion. Consequently, the 
following question can be immediately raised. Since the torsion or twist of the 
blade is assumed to occur about the elastic axis of the blade which can have a 
different orientation as the blade flaps and lags, the spring representing the 
torsional stiffness of the blade should also change its orientation! The answer 
to this question is negative. Since the model is intended to represent a hingeless 
cantilevered blade for which the slope of the elastic axis at the root is always 
zero, irrespective of blade orientation. In our model, the torsional stiffness of 
the blade is represented by an equivalent torsional spring at the root. Hence, its 

orientation does not change with the blade motion. The spring stiffness are Kn for 

B 

flaping, K r for lead-lag, for torsion and Ka for the pitch link flexibility. 

From the physics of the problem it is clear that these springs are all torsional 
type springs. 

In deriving these equations the sign convention is important. All counter 
clock-wise rotations and moments are taken as positive. The restoring moment 
in any torsional spring, due to a positive displacement (i.e., rotation), is 
clock-wise and hence negative. 

Another important ingredient in this derivation is the coordinate system. Let 
x,y,z be an orthogonal triad attached to the undeformed blade with zero pitch angle. 
The X-axis is along the elastic axis of the blade as shown in Fig. A2 . The X’, Y', 
Z' system represents another orthogonal triad, rotated through an angle 0 about 
X-axis in the counter clock-wise direction. The angle 0 represents the collective 
pitch of the blade shown in Fig. A2 . The root springs representing the blade flexi- 
bility are oriented as follows: Kg along Y' axis, Kj- along Z' axis and K^ , K^- 

B B B C 

along X' axis. 

Using this information the expressions for the moments are derived next. 

To derive the expressions for the restoring moments of the springs, due to 
blade motion, the total angular displacement of the blade has to be decomposed 
into components along the directions of the spring restrained hinges. The rotation- 
al (angular) displacement components are then multiplied by the corresponding spring 
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stiffness to yield the appropriate moments about the hinges. Finally, the various 
moment vectors are expressed in terms of the component acting along the undeformed 
blade axes X, Y and Z respectively. These expressions are then compared with the 
results obtained in previous studies . 

To arrive at the deformed blade orientation from the undeformed position a 
specific sequence of rotations are followed, namely flap, lag and torsion. The unit 
vectors in the directions of the undeformed and deformed blade coordinate systems 
are related by an Euler angle transformation which is derived below. 

The flap rotation is assumed to take place first, thus the X,Y,Z axes system 
is rotated through an angle (-3), in the clockwise direction, about Y-axis repre- 
senting the flap hinge of the blade shown in Fig. A3. Components of the angular 
displacement (-3) along X', Y', Z' axes can be obtained using the transformation 



0 0 

cos0 sin0 

-sin0 cos0 



(A . 1) 


where e , e ,, e , and e , e , e are unit vectors along 
x y z x y ’ z 

respectively . 

Components of (-3) along X', Y', Z' axes are along 


X’, 


Y\ Z' and X, Y, Z axes 


X' - 0 
Y’ - -3cos0 

Z' - 3sin0 (A. 2) 


Due to this rotation (-3) , the coordinate system is rotated to a new position 
X^, Y^, Z^ shown in Fig. A3. The transformation of unit vectors between the two 
systems is given by 


k] 

1 

cos3 

0 

-sin3 

P ! 


0 

1 

0 

k ) 


sin3 

0 

co s3 



(A. 3) 


The lag rotation is assumed to take place next, the system X^, Y^, Z^ 
is rotated about Z x axis through angle (?) , in the counter clock-wise direction, 
representing lead-lag motion of the blade, shown in Fig. A. 4. Equation (A. 3) is 
used to obtain the components of the angular displacement '?' along X,Y,Z axes. 
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These components are 


along X - - £sing 
Y - 0 

Z - £cosg (A. 4) 

These components along X,Y,Z are transformed into components along the dir- 
ection of the spring axes x'»Y'»Z' using equation (A.l). These components 
are: 


along x' - -£sin£ 

Y* - CcosgsinQ 

Z 1 - CcosgcosO (A. 5) 


Next the X^, Y^, Z^ axes system is rotated to a new position X 2 , Y£, Z^ 
shown in Fig. A4 . The relation between the unit vectors along X 2 , Y 2 , Z 2 
and Xj, , Y-^, Z-^ systems is given by 


( §xl 


cos t, 

-sinC 

0 

V 


sinC 

cos? 

0 

w 


0 

0 

1 



(A. 6) 


Combining equations (A. 6) and (A. 3), the transformation of unit vectors between 
X^ , Y 2 > 7^ anc * X, Y, Z systems; given below is obtained. 


e 

x 


e 

y 


e 

z 


cosgcosC 

-cosgsini; 

-sing 

( 6x2 ) 

sin? 

cosC 

0 

1 ty2 1 

singcosC 

-singsinC 

cosg 

( e z2 ) 


(A. 7) 


Finally the torsional rotation is assumed to take place. To represent torsion, 
the blade is rotated through an angle (<j>) in the counter clock-wise direction 
along the X 2 axis, as shown in Fig. A5. The components of the angular dis- 
placement ((f)) along X, Y, Z system are obtained using equation (A. 7) . These 
components are respectively: 


along X - <j>cos£cos£ 
Y - <j>sin£ 

Z - (fisingcosC 
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(A. 8) 



These components are again transformed in the X', Y', Z' directions using 
equation (A.l). The corresponding components are 

along X' - (JicosgcosC 

Y' - <j)sinZicos8 + <j)sin0sinBcost; 

Z' - -<}>sin?sin0 + <j>cos0sinBcosC (A. 9) 

The total twist of the various springs oriented along X', Y', Z' axes is 
obtained by adding equations (A. 2), (A. 5) and (A. 9). The various components 
can be identified as: 

along X' - -£sin0 + <J>cos0cos£ 

Y' - -geos© + £cos0sin0 + c|>sinZ;cos0 

+ <J>sin0cos£sin0 

Z' - 0sin0 + £cos0cos0 - <J>sin£sin0 

+ 4>sin3cos£cos0 (A. 10) 


In the model shown in Fig. Al, the springs representing the torsional 

stiffness of the blade and the pitch link stiffness are in series and are 

along X' axis. These two springs can be combined and can be represented by 

an equivalent spring of stiffness K, = K , - K. / (K. + K . ) , when K, is very 

9 9 B q>c ‘Pb ‘Pc <Pq 

large, thenK, becomes K, . The restoring moments in the springs, due to 

<p cpg 

the flap, lag and torsional rotation of the blade can be written as 


M^, = -K^ (-£sin0 + (JjcosBcosC) 

M ' = -K„ (-0COS0 + £cos0sin0 + (f>sin£cos0 + 4>sinBcos£sin0) 

y p b 

K , = -K_ (Bsin0 + £cos0cos0 - <{)sin^sin8 + 4>sinBcos£cos8) (A. 11) 

z C B 

Transformation of these moments alone the undeformed blade axes (X,Y,Z), yields 
the torsional, flapping and lead-lag restoring moments, associated with blade 
motion. These are: 
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Torsional moment 


M = - K, (-£sin8 + (jicosgcos?) 
x <P 


Flapping moment 

M. = -Ka cos0 (-8cos0 + Ccos8sin0 + <5>sin^cos0 + <t>sin8cosCsiu8) 
y P B 

+ K_ sin0 (8sin0 + £cos8cos0 - <j>sinCsin0 + <j>sin8cos£cos0) 

Lead- lag moment 

M = -Kg sin8 (-0COS0 + £cos8sin0 + <j>sin£cos0 + <|>sin8cos£sin0) 
z P B 

-Kr cos0 (0sin0 + ^cosScosQ - <t>sin^sio.0 + <t>sin.6cosCc.os9) 

B (A. 121 

Assuming 8 and Z, are small angles and making the approximation sin8 - 8, 
sin^ = £, cos8 - 1 and cos£ - 1, equation (A. 12) can be simplified to yield 


M 

y 


M 

z 


-K^C-58 + 4>) 

-Kg cos0 (-8cos0 + £sin0 + <j>£cos0 + (J>8sin0) 

+ K r sin8 (Bsin8 + £cosB - 4>£sinB + <)>BcosB) 
^B 

-Kq sinO (-0cos0 + CsinO + <Hcos0 + <|>Bsin0) 

Pb 

-K r cqs0 (8sin0 + £cos0 - <J>£sin0 + cj>Bcos0) 


(A. 13) 


Rearranging the terms, 

M = -K, (-C8 + <J>) 
x <j> 

M = (8 - <j>£) (Kg cos 2 0 + Kj- sin 2 0) + (£ + <J>8) (K^ - Kg ) sin0cos0 

y B B 3 B 

M = -(C + 4>B) (Kr sin 2 0 + K r cos 2 0) - (8 - 4>C) (Kr - Kg ) sinOcosB 
z P B 5 B ^ 3 B 

(A. 14 

Equation (A. 14) represents the torsional, flap and lag restoring moments of 
the springs due to the flap, lag and torsional motion of the blade. 

These expressions for the restoring moments acting on the spring re- 
strained model of the blade are compared with results published in previous 
studies. For completeness these comparisons are carried out for both spring 
restrained models of a hingeless blade and fully elastic representations of 
th£ hingeless blade. 
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Ormiston and Hodges [Ref. 13] presented moment expressions for a spring 
restrained model of a rotor blade representing a hingeless blade. They have 
also taken hub flexibility into account. In order to compare the results 
obtained in Ref. 13 with the expression obtained in this report, the hub stiff- 
ness is allowed to become infinite, thus only the simplified expressions are 
compared. Recall that the torsional degree of freedom was not considered 
in Ref. 13. The expressions [Ref. 13] for the restoring moments are: 

Flapping moment (Eq. 44 of Ref. 13) 

M = -3(Ko cos 2 0 + K r sin 2 0) - £(Kr - K Q )sin0cos0 (a) 

y » b S 

Lead-lag moment (Eq. 45 of Ref. 13) 

M = -£(K fl sin 2 0 + K r cos 2 0) - 3(Kr - Ko ) sin0cos0 (b) 

z ^B S S B 

Setting <f> = 0 in Eqs. (A. 14) ' yields 
M x = K <f> R » 

M = B(K q cos 2 0 + K r sin 2 0) + £(K_ - Ko ) sin0cos0 

y PB P B 

M = -?(K r sin 2 0 + K r cos 2 0) - 3(K r - Ko ) sin0cos0 (A. 15) 

z P B '’B P B 

By comparing equations (a) and (b) with equation (A. 15), we find there is a 

negative sign in the expression for M . The reason for this discrepancy is 

that the authors of Ref. 13 used different sign convention. In their 

case. M clockwise is positive and M counter clock -wise is also positive. 

The authors of Ref. 13 have used this sign convention because their objective 

was to derive a blade model which was a simple analog to the flexible blade 

equations (hingeless) which were derived by Houbolt and Brooks [Ref. 15] who 

have employed the same sign convention for M and M . The expression of the 

y z 

present study are also compared with the results of an elastic blade model. 

The moment expressions obtained by Rosen and Friedmann [Ref. 16] for an 
elastic blade are (Eqs. (15) of Ref. 16) 


M 

x 

M 

y 


G J ( c[) + v w ) 

,x ,xx ,x 


= -(EI 2 - EI 3 ) sin9cos0(v xx + <J>w^ xx ) 
-(EI 2 Sin 2 0 + EI 3 cos 2 0)(w xx - <f>v xx ) 
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M = (El - El- ) sin0cos0(w - <j)v ) 

Z Z j ) XX y XX 


+ (El- cos^G + El- sin^0)(v + <]>w ) (c) 

/ J j XX y XX 

where v, w, 4> represent the elastic lead-lag displacement, flap displacement 

and twist at any section of the blade, and the quantities EI^, EI^ , GJ 

represent lead-lag, flap and torsional stiffness of the blade. 

The expressions obtained in Ref. 16 represent the applied moments, thus 

to obtain the restoring moments Equation (c) have to be multiplied by (-1) 

and the resulting expressions are compared with Eqs. (A. 14). The expressions 

for M and M are the same, however there is a sign difference in the ex- 
y z 

pression for torsional moment. The expression derived in this study contains 
a term (4>-£f3) while the corresponding term in Ref. 16 is (4> + v w ). 

J X | XX y X 

This discrepancy is due to a different sequence of rotation followed by Rosen 
and Friedmann [Ref. 16]. The sequence of rotations adopted in Ref. 16 was 
lead-lag, flap and torsion. When a sequence of flap, lead-lag and torsion, 
such as employed in the present study is used, the results of Ref. 16 are in 
agreement with the results obtained in this study. Equations (c) , representing 
a hingeless blade, with Eqs. (A. 14) corresponding to the rigid, offset hinged 
spring restrained blade model one can identify a number equivalence relations, 
which provide some physical insight. These equivalence relations are given 
below. For torsion one has 


GJ(<f) ,x - V ,xx w .x ) = V* - «> 


(A. 16) 


where 

K, = K , K, / (K, + K, ) 

4- 4> b <f> c * B 4> c 

If the pitch link flexibility is very large, i.e. -*■ 00 then becomes 
K ' 




We shall consider only K. in Eq . (A. 16) because Rosen and Friedmann have 


no? considered the pitch link fSexibility. Thus equation (A. 16) should be 
rewritten as 


GJ( ^,X - v ,xx w ,x ) = K «0 b ( * - 
Comparison of the various other terms yields 

<E1 2 - EI 3 )( ”,» + ' «r - + »> 

B B 

EI- (w - <|>v ) = K (B - 4>c) 

2 ,xx T ,xx r 

B 


(A. 17) 

(A. 18) 
(A. 19) 
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(A. 20) 


E V w ,xx - *\x*> = K e B (e " 


(El - El )(w - 4>v ) = (K. - Kg )(B - <frO 

Z J ,XX )XX Sg Mg 


EI 2 (v ,xx + ^.xx 1 - K ? r « + ♦« 


EI 3 (v ,xx + <f,w ,xx ) = \ + ♦« 


(A. 21) 
(A. 22) 
(A. 23) 


Setting c{) = 0 in equations (A. 19), (A. 20), (A. 22) and (A, 23) yields 


El, w = Kr B 
2 , xx G R 


El, w = K ft g 
3 ,xx P B 


EI 2 \xx = V 


El, v = K„ ? 
3 ,XX Pr 


(A. 24) 


Examination of these expressions reveals that Kj- , Kg represent the lead- 

B ' B 

lag and flap stiffness EI 2 , EI^ respectively and g,£ correspond to the 
curvatures in the elastic blade analysis. Comparison of the appropriate terms 
in equation (A. 17) reveals 


GJ = K, 


B 


<t> 

,x 


= <f> 


w 


,X 


= e 


(A. 25) 


This comparison indicates that both w and w can be identified as 
v ,xx ,x 

quantities corresponding to the flapping angle 3 in the rigid spring 

restrained blade model. This statement requires further clarification. 

The comparison outlined above indicates that when dealing with moment 

terms, such as El w ,3 has the role of the curvature. However when 
’ ,xx’ 
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examining the role of 3 in the various transformations relating the unit 
vectors of the deformed and undeformed states of the blade the role of 3 
corresponds to w in the flexible blade equations. Furthermore it can be 

i X 

seen that <p in an elastic blade corresponds <j) in a rigid blade model, this 
> x 

is due to the fact that the torsional moment in an elastic blade is defined 

as GJij) whereas in the spring restrained blade model the torsional elastic 
> x 

restoring moment is K^<}>. 

A. 3 Moment Equations for an Articulated Blade 

A typical articulated blade, which has no root springs in flap and lag 

will not experience elastic restoring moments in about these hinges. Therefore 

the model proposed for this case, shown in Fig. A6, involves only a torsional 

spring and no springs for flap and lag. In the model shown in Fig. A6, two 

torsional springs are provided. The spring K, represents the blade torsional 

<Pb 

stiffness and Kj represents the control link stiffness. These are aligned 
along the undeformed blade axis x. Since the blade can perform flap and lag 
motions during operation, the spring Ka can orient itself along any direction 
in space provided the pitch link is inboard of flap and lag hinges. The spring 
remains along the uhdeformed x-axis. If the pitch link is outboard of 
leading and flap hinges, then both and change their orientation as the 
blade undergoes flap and lag motion. Therefore one needs to consider two 
cases: (a) pitch link inboard of flap and lag hinges and (b) pitch link out- 

board of flap and lag hinges. 

Case (a) Pitch Link Inboard of Flap and Lag Hinges 

Consider Fig. A. 6 numbers 1, 2, 3 refer to three nodes. At node 2, the 

two springs and Ka are connected. The springs are oriented along the 

B C 

undeformed elastic axis (X-axis of the coordinate system) . In this case, since 

the pitch link is inboard of flap and lag hinges, only Ka takes on different 

orientation as the blade undergoes flap and lead-lag deformation. 

Assume that the blade has undergone a flap deflection (-3) and a lag 

deflection (£) respectively, as shown in Fig. A7 . The new position of the 

elastic axis is and the deformed coordinate system is X£, Y ^ , Z The 

spring associated with blade torsional stiffness, Ka , is now oriented along 

B 

the X 2 axis and the spring associated with pitch link stiffness, K^ is oriented 
along the undeformed axis, i.e. X-axis. For this case one may write a stiffness 
matrix relating moments to angular rotations as indicated in Eq . (A. 26). 
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Ka cosBcosCsinC 

sin 2 C 

B 



Symmetric 



2 

Ka cosBcos CsinB 
V B 

Ka sinBcosCsinC 
B 

Ka sin 2 Bcos 2 £ 
b 





-Ka cos 2 Bcos 2 £ 

-Ka cosBcosCsinC 
y B 

2 

-K^ cosBcos CsinB 
^B 

Ka +Ka cos 2 Bcos 2 C 




-Ka cosBcosCsin? 
^B 

-K* sin 2 ^ 

V B 

-Ka sinBcos^sinC 

Ka cosBcosCsinC 

Ka sin 2 C 
B 


il 

2 

-Ka cosBcos CsinB 
B 

-Ka sinBcosCsinC 
B 

-Ka sin 2 Bcos 2 C 
B 

2 

Ka cosBcos CsinB 
B 

Ka sinBcostsinC 

2 2 

Ka sin Bcos t 
B 


0 

0 

0 

-Ka 

‘•‘c 

0 

0 

K0 

C 

0 

0 

0 

0 

0 

0 

0 0 

0 

0 

0 

0 

0 

0 

0 0 


(A. 26) 



This stiffness matrix relates the generalized displacements (rotations along X, 
Y,Z directions) and the moments at nodes 1,2,3. In equation (A. 26), ^ re- 
presents the moments along X,Y,Z axes at node i and 1 represents the ro- 

x , y , z 

tations along X,Y,Z axes at node i. 

The equivalent stiffnesses of the system in torsion are obtained from the 
solution of this matrix equation. The conditions under which this equation is 
solved are 


<T = C 


<(> 3 = 0 
z 


and 


M 


2 2 

M = M = 0 
y z 


Node 3 is fixed 


No external moment at node 2 


When trying to solve the matrix equation under these conditions the individual 
equations for M^, M^, become redundant. This redundancy is due to the rigid 
body degrees of freedom being included in the model. The physical reason for this 
redundancy is due to the fact that the moment along the Y and Z axes cannot be 
resisted by this spring model. In order to overcome this difficulty, it is as- 
sumed that the flap and lead-lag angles are very small while the spring K. is 
oriented along the undeformed elastic axis, i.e., X-axis. The springs Ka and 
K.^ are in series and are oriented along X-axis and can be combined. For this 
model, the relation between the torsional moment and the twist <J) of the blade 
becomes 

M = [K. K. /(K. + K. ) 34 > 


= K . 


(A. 27) 


Case (b) Pitch Link Outboard of Flap-Lag Hinges 

When the pitch link is outboard of the flap-lag hinges both the springs 
and are always oriented in the same direction, from which the following 

refation between torsional moment and elastic twist is obtained. 


M 

x 



K, 

<p r K 


+ V'* 


(A. 28) 


Therefore for an articulated blade, notwithstanding whether the flap and lag 
hinges are inboard or outboard of the pitch link, the relation between the 
torsional moment and the elastic twist is given by 


M 


[K* K a /(K, + 


4> b 


B 


k m 


= v 


(A. 29) 
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Combining this result with the result for a hingeless blade, discussed 
previously, the general elastic restoring moment expressions become; 


Torsion: = - K^ (4> - £3) - K^ 2 <}> 

Flap; M = (3 - <f>?)(K R cos 2 8 + K c sin 2 0) + (({>3 + ?)(K C -Kg )sin6cos0 

y P B B B P B 


Lead-lag: M = -(£ + <t>3) (Kr sin 2 0 + Kr cos 2 0) - (3 - 4>£) (Kr -Kg )sin0cos0 

z P B B ^B B 

(A. 30) 


For a hingeless blade 


K$ = 0 and 


V*. 


for an articulated blade 

\ " °> k 3 b = °* % = 0 
and K 4 , c K <j)B 

K * 2 ' \*\ 


A. 4 Elastic Restoring Moments on a Rigid Blade with Root Springs and Hub 
Flexibility 

In this section, the elastic restoring moments due to the root springs 
and due to the hub flexibility are derived. The model for this system based 
on equivalent springs is illustrated in Fig. A8 . There are two sets of ortho- 
gonal springs, springs with stiffness Kg , K r and Kj which represent the 

P B 

blade stiffness in flap, lead-lag and torsion, respectively. The spring 

constant K^ represents the stiffness of the pitch link system and Kj- and 
C H 

Kr are the stiffnesses of the hub in flap and lead-lag directions. The hub 
P H 

is assumed to be torsionally rigid. The X,Y,Z axis system represents the 
undeformed blade coordinates with the X-axis oriented along the elastic 
axis. The system X', Y', Z' corresponds to the X, Y, Z axis system after 
a counter clock-wise rotation about the X-axis by an angle 9. The angle 
0 represents the collective pitch of the blade. 
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The orientation of the various springs is as follows: K is along 

Z-axis, K 0 is along Y-axis, K c is along Y'-axis, K is along Z'--axis, 

6 H 6 B C B 

K, and K, are along X-axis (see Fig. A8) . The numbers 1, 2, 3 in Fig. 

A8 refer to the three nodes. The blade spring system is attached to the 

hub spring system at node 2. 

The relation between the moments and the angular displacements for this 
equivalent spring system describing the combined blade and hub is given by 
the matrix equation (A. 31), on the next page. 

In equation (A. 31) M*, M^, refer to the elastic moments at node i 
in the X,Y,Z directions, respectively and <J>^ refer to the angular 

displacements at node i in the X,Y,Z directions, respectively. The various 
moments are obtained in terms of the angular displacements by solving this 
matrix equation subject to the following conditions: 

When the hub is fixed, the angular displacements at node 3 are zero and then 


3 3 3 

b = d> = d) = 0 
r x ^ y r z 


(A. 32) 


Since there is no external moment at node 2, 


2 2 2 

M = M = M = 0 
x y z 


(A. 33) 


Solving the matrix equation, Eq. (A. 31) subject to the conditions given in 

2 2 2 

equations (A. 32) and (A. 33) and recognizing that = 0 one 

obtains 


2 2 1 1 
-(K cos 9+K c sin 0) (p (K. - K ) sin0cos0<j> 

h P B Z 3 B C B y 


-(K + K cos 2 0 + K_ sin 2 8). (j) 2 - (K D - K ) sin0cos0<J> 2 

S ' 6 B 2 6 B C B y 


(A. 34) 
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(A. 31) 



•(Kg - Kr )sin0cos0 <J>^ - (Kg cos 2 8 + Kr sin 2 0) (J) ^ = 
B B 2 B B^ 


- 4> (Kg - Kr ) sinScosS- (Kg + Kg cos 2 0 + K r sin 2 0)<fi 2 

Z Pg ^g P H Pg <*g y 


- \ +x ■ -«♦, + V *1 


(A. 35) 
(A. 36) 


2 2 2 ' 

From Eqs. (A. 34), (A. 35) and (A. 36) <}> and <j>^ can be expressed in terms 

- ,1 A ,1 

of <p , <p , cp , thus 
x y z 


>1 • V^B + V 


(A. 37) 


- 


- 

- 


( A 

*1 


Kg +Kg cos 2 0 
e H P B 

-(Kg -Kr ) 
P B ^B 


(K r cos 2 0 + Kg sin 2 0)4>^ | 

S P B 2 1 



+Kr sin 2 0 

S 

sin 0cos0 


+(Kg -Kg )sin0cos0$J I 



-% - V 

Kr +Kr COS 2 0 
H B 

> 

l ^ 

■ (Kg -K r )sin0cos0<}> 1 
P B 2 

- - 


sin0cos0 

+Kg sin 2 0 

B J 

1 

+ (Kg cos 2 0+K r sin 2 0)4>^ 

, ^B ? B y 


1 

A 


where 


(A. 38) 

A = (Kj- +K^cos%+Kg sin%)(Kg +Kg cos^0+K^ sin 2 0)-(Kg -K^ ) 2 sin 2 0cos 2 0 
Hd b hb b bb 

= K r Kg +K r (Kg cos 2 0+Kr sin 2 8)+Kg (Kg sin 2 0+K r cos 2 0)+K r Kg (A. 39) 

P H S P B T5 P H p b < *B ^B P B 

From Equation (A. 38) 

,2 _ 1 


4> - T" {[Kg (Kr cos 2 0+Kg sin 8)+Kr Kg ]4>^ + Kg (Kg -Kr )sin0cos0<j>"'"} 

Z A P H P B P B 2 P H B B v 


y 

(A. 40) 
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and 


4>y = |-{K ? ( K B B “K^)sin0cos0^+ [K^ (Kg^cos 2 6+K ? sin 2 0)+Kg K^] <f>J } (A. 41) 


2 2 2 111 

Substituting for <j> , < j> , <J> in terms of $ , (|> , 4> in equation (A. 31), and writing 
y J 2 , x , yzx 

the expressions for M , M , M , one obtains the torsional moment 


K<}, 

H 1 - - B - C - / 

x Ka +Ka V x 


(A. 42) 


The flap moment can be written as 

,1 


M J = f {1 ^ H Ke H (K 6 B COs2e+ S 8in : 6) + W H n 


- -r- (K r -Ko ) K r Ko cosBsinB 

A ? B e B C H P H 


(A. 43) 


The lead-lag moment can be expressed as 


1 t 2 2 

M = — {K r K q (K, cos 0+K ft sin 0) + K c Kr Kr 

z A C H 3 H V C b 3 b ? h b B 


- -jf- (K r -Kr ) Kr Kr sin0cos0 
A ^B P B P H 


(A. 44) 


In the last three equations <f)^, tf)^, (f)^ represents angular displacements 
given at node 1 along X,Y,Z axis. These angular displacements can be related 
to the flap, lead-lag and torsional rotations of the blade. The flap angle 
(-3) is about the undeformed Y-axis, the lead-lag angle ( cj)) is about the Z axis 
after it has undergone flapping deformation and the torsion (£) is along X-axis 
after the flap and lag deformations have taken place. Resolving these angular 
displacements along X,Y,Z axis, using equations (A. 4) and (A. 8) together with 
the fact that (-3) is along Y-axis, the components of the angular displacement 
(or rotation) about the X,Y and Z axis are 


<!>^ = -£sin3 + <J)cos3cos£ 

4>y = -3 + <j>sin£ 

<j>^ = £cos3 + $sin3cos£ (A, 45) 
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Introducing the small angle assumptions for g, £ and , equations (A. 45) become 

= - £3 + 4> 


** = - 3 + 


= ? + 4>B 


(A. 46) 


Substituting these in equations (A. 42), (A. 43), and (A. 44) , the elastic moments 
are obtained 


M L = (K r K R (K r cos 2 0'+ Ko sin 2 0) + K r Ko K r } 

z A ? H 0 H C H P B S 


(3 " ^ C) (K r - K fi ) K, K, 


?B - V 


sin0cos9 


(A. 47) 


M = 


M 


- { Kc Kfi __ ( Kr cos 2 0 + R r sin 2 0)+K R K r K R } 

am) 

A 


. >.K r Ko 1K R cos o t k_ sin o;-t-K R *. r ^ R 

A P H P B ^B P B ^B P H 


. (K r - K r ) K r K r sin0cos0 
A ^B P B P H 


B r C 
K. +K , 

+ B V 


(<f> - 


(A. 48) 


(A. 49) 


Equations (A. 47), (A. 48) , and (A. 49) can be written in an alternative form, i.e. 
in terms of equivalent stiffness in flap, lag and torsion respectively. 


Defining 


K, 


K„ 


B C 

V^c 


H K g 

H B 
Kg +Kg 
H B 


H B 
K ? +K ? 

H B 
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equations (A. 47), (A. 48), (A. 49) with the superscript 1 deleted, can be 
written as 


M 

z 


M 

y 


M 

z 


{K - R(K - K ft ) sin 2 0} + 

A' ? C P 

- {K ft + R(K - K ft ) sin 2 0} 

A , e ? e 

(<f) - ?B) 


IMP 

A ' 

(C-H»B) 

A' 


R (K 
- R (K^ 


V 

■ y 


sinGcosG 


sin0cos0 


(A. 50) 
(A. 51) 
(A. 52) 


where 

1/Kg - 1/K ? 


(1/Ko + 1/Kg ) - (1/K C + 1/Kj- ) 

P B P H 

and A' = 1 + R (1 - R) sin 2 0 (K - K g ) 2 /K & K ? 


When the hub stiffness are very large R becomes unity and the moment equations 
reduce to those obtained previously. When the blade stiffness is very large, R 
reduces to zero. These moment expressions have to be multiplied by (-1) to 
obtain the restoring moments of the springs. 

Peters [Ref. 14] has given the moment expressions for a spring model 
including hub flexibility. The torsional degree of freedom was not consid- 
ered in Ref. 14, thus when substituting <J> = 0, equations (A. 50) - (A. 51) 
reduce to the equations given by Peters. 


Combining the results for an articulated blade and for an elastic blade, 
general expressions for moment are given below. 

The restoring moment in flap is 


M = {k r + R(K - K c ) sin 2 0}+ 


y a' e 

The restoring lead-lag moment is 


A' 


R(K - K d ) sin0cos0 

4 P 


(A. 53) 


\ = - ( ^ g) (K r - R (K r - K r ) sin 2 0} - <t>g) R(K ? -Kg) sin0cos0 


A' 

The restoring torsional moment is 


A’ 


(A. 54) 


M X = * \ ^ ♦ 


(A. 55) 
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For an articulated blade 


K 4> % 

K 6 - K - K* - 0 ..X K - ^ - 

1 2 *c *8 


For a hingeless blade K, =0 


; \\ 

6 " v\ 


K - ^ 

' " \ *\ 


1/K. - 1/Kr 

21 J 

(1/K r + 1/K ft ) - (1/K + 1/K ) 

B B ? B ? H 


and A' = 1 + R (1 - R) sin 2 0 (K - K Q ) 2 /K r K 0 

Q> P C P 
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Figure Al. Equivalent Spring Restrained Blade Model 




z 



Figure A8. Equivalent Spring Restrained Blade Model With Hub 
Flexibility 
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